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… Klassieke mechanica — 
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De klassieke mechanica, zoals hier behandeld, beoogt een voort- 
zetting en aanvulling te bieden van de voorkandidaatscolleges klassieke 
mechanica. De onderwerpen die ter sprake zullen komen zijn: principe 
van d'Alembert, de vergelijkingen van Lagrange, het principe van 
Hamilton, de vergelijkingen van Hamilton, de theorie van kleine 
trillingen en de viriaalstelling. 

Om te beginnen herinner ik nog even aan de wetten van Newton (1686), 
die het uitgangspunt vormen van de klassieke mechanica, Deze grondwetten 
berusten uiteraard op ervaringsfeiten, al is het niet eenvoudig con- 
crete experimenten aan te geven, waarop ze gebaseerd zijn. Ook hier ís 


iet feitelijk weer de goede overeeristemming van de gevolgtrekkingen met 


‚de ervaring, die een rechtvaardiging vormt. In zijn Foundations of 


Physics waarschuwt Frank vooral tegen de mening, dat de wetten van Newtcn 
min of meer vanzelfsprekend zouden zijn. Hij zegt, _ meer men ze als 
evident beschouwt, des te minder heeft men ze begrepen. Dat de wende 
ten niet zo voor de hand liggen, als het nu wel eens ons voorkomt, blijkt 
zak duidelijk uit de historische ontwikkeling van de mechanica van 


Aristoteles tot Newton (vgl. hiervoor het boek "Val en Worp" van 


Dr. E.J. Dijksterhuis). Aristoteles beschouwde de vrije val als de natuur- 
lijke beweging, alle andere bewegingen noemde hij gedwongen. Voor Galilei 


‘meende men ook veelal, dat voor het onderhouden van een constante snel- 


heid al een oorzaak (kracht) nodig was, niet alleen voor het geven van 


"een versnelling. 





Newton's grondstellingen worden ook niet meer critiekloos aanvaard. 
Men denke bv. aan de relativiteitstheorie van Einstein, waarin de be= 
grippen absolute ruimte en tijd van Newton verdwenen zijn. We zullen 


ons verder beperken tot enkele opmerkingen over de wetten afzonderlijk. 


1. De traagheidswet (feitelijk reeds afkomstig van Galilei). 
Formuleren we deze als volgt: 


Een massapunt (onder een stoffelijk of massapunt verstaan we een 
geometrisch punt met massa), waarop geen uitwendige invloeden werk= 
zaam zijn, blijft ín rust of in eenparige rechtlijnige beweging, 


dan zegt deze al heel weinig, omdat er altijd wel een coördinatenstel- 
sel is aan te wijzen t.o. waarvan het massapunt in rust is of eenparig 
beweegt. Een betere formulering is al: 


Als voor een bepaald massapunt, dat niet onderworpen is aan uitwendige 
invloeden, geldt dat het t.o.v. een bepaald coördínatenstelsel in 

rust of in eenparige beweging ís, dan zal elk ander massapunt, dat 
vrij is van dergelijke invloeden (d.w.z. dat ver genoeg van andere 
massapunten verwijderd is) zich t.o.v. ditzelfde systeem eenparig 
rechtlijnig voortbewegen of ín rust bevinden. 


Eigenlijk postuleren we dus dat een dergelijk stelsel bestaat. Zo'n 
coördinatensysteen heet een Snerttanlstelsel. Het is duidelijk, dat 
ieder stelsel dat t.o.v. een inertiaalstelsel eenparig rechtlijnig 
voortbeweegt zelf ook een inertiaalstelsel is. Belangrijk is, dat we 
voor vele doeleinden een vast met de aarde verbonden stelsel als iner= 
tiaalstelsel kunnen beschouwen. Niet altijd is dit juist, zoals bv. 
blijkt uit de slingerproef van Foucault. Een betere keuze is dan een 
stelsel met de oorsprong in de zon en waarvan de assen vaste richtingen 
hebben t.o.v, de sterren. n 


In onze formulering dreigt verder nog gevaar, dat het traagheids- 


beginsel tot een tautologie wordt, als we nl. van "vrij van uitwendige 











invloeden" eerst zullen spreken als:het beschouwde massapunt eenparig 
beweegt of in rust is. Iets dergelijks kan optreden bij tijdmetingen: 
zouden we een tijdinterval definiëren als evenredig met de hoek waarover 
de aarde gedurende dat tijdinterval draait, dan wordt de bewering dat 

de aarde eenparig om zijn as wentelt tot een tautologie en wordt de 
discussie over de afwijking van de aardrotatie van een eenparige rotatie 
paradoxaal. Dit laatste probleem wordt slechts zinvol omdat men later 


de tijdmeting anders is gaan definiëren. 


2, Voor een verandering van de snelheid t.o.v. een inertiaalsysteem ís 
een oorzaak nodig, die we kracht noemen; de samenhang tussen kracht en 
versnelling zal zijn: 


md 
R-ma (of beter = a) 


waarin de evenredigheidsconstante m de trage massa voorstelt. Hier hebben 
we de moeilijkheid dat KX noch m nauwkeurig gedefinieerd is. Newton defi- 
nieerde m als volume maal dichtheid, maar dit leidt kennelijk tot een 
vicieuze cirkel, immers de dichtheid is weer massa gedeeld door volume. 
Men zou bv. krachten onafhankelijk met benap van uitrekkingen aan veren 
kunnen definiëren. Het blijkt dan, dat we krachten als vectoren kunnen 
samenstellen (parellelogram van krachten). Een kracht is een vector, díe 
dezelfde richting heeft als de door hem veroorzaakte versnelling. Verder 
blijkt dan voor elk lichaam de verhouding van kracht en versnelling een 
voor dat lichaam karakteristieke constante te zijn; de zg. trage massa, 
die bij homogene substanties evenredig is met het volume. De eenheid van 
massa is de kilogram. De kracht die hieraan een versnelling geeft van 
1 mfsec* heet de Newton (= 10° dyne) . | 

Uit nauwkeurige metingen is gebleken dat op een bepaalde plaats op 


aarde lichamen met dezelfde trage massa door de aarde even sterk worden 





aangetrokken of algemener, dat het gewicht van een lichaam (d.i. de 
aantrekkingskracht door de aarde er op uitgeoefend) evenredig is met zijn 
trage massa. Voor elke plaats op aarde is er dus een karakteristieke 
constante, de versnelling van de zwaartekracht ter plaatse. Men drukt 

dit soms uit door te zeggen: zware massa = trage massa. Einstein heeft 
als eerste op de merkwaardigheid van dit feit gewezen. De zwaartekracht 
is de enige kracht, die evenredig is met de massa van het lichaam waarop 
hij werkt (de electrische kracht bv. ís er geheel onafhankelijk van). 

Het gevolg is dat we niet kunnen uitmaken of we in een zwaarteveld zitten 
of wel in een stelsel dat een zekere versnelling t.o.v. een inertiaal- 
stelsel heeft. Dit zgn. aequivalentieprincipe is het uitgangspunt geworden 


voor de algemene relativiteitstheorie, 


35. Als derde wet van Newton beschouwt men gewoonlijk de wet: actie = 
reactie. D.w.z. de krachten die twee massapunten (bv. tengevolge van 
hun gravitatie of van hun electrische ladingen) op elkaar uitoefenen, 
zijn gelijk en tegengesteld, dus x, + K= O. Het ís daarbij niet noode= 
zakelijk, dat deze krachten langs de verbindingslijn der massapunten 
zentent eid, (Bv. de krachten die bewegende ladingen op elkaar uitoefe- 


nen.) Ook deze wet moet gezien worden als ervaringsfeit. 


We gaan nu verder uit van de geldigheid, t.o.v. een inertíiaalstelsel, 
van de wet rd = mi, waarbij we in vrijwel alle gevallen, tenzij uit- 
drukkelijk vermeld een aan de aarde vast verbonden stelsel als zo'n 
inertiaalstelsel zullen beschouwen. 

Als nu de krachten als functies van de coórdinaten en de tijd 
(eventueel ook van de snelheden) gegeven zijn - de krachten aan andere 


gebieden der fysica ontleend =- vinden we de beweging door de 3 differen- 














tiaalvergelijkingen van de 2° orde 

ee mi(t) = K(x,y,z,t) 

u my(t) = K(x,yzot) 

ä mz (t) = K‚(x,y,z,t) 

ie te integreren. De oplossing is volkomen bepaald als op een gegeven 
tijdstip to plaats en snelheid bekend zijn. 

pie We moeten wel bedenken dat de 2° wet van Newton zonder meer slechts 

) geldt voor een inertiaalstelsel. Voor een stelsel Ss, 0902) dat 

Ed 

4 eenparig beweegt t.o.v. een inertiaalstelsel S(x,y,z) is 

df 

Ad a, =e 

d 1 

8 zodat weer geldt FE = ma, . Ook S, is immers een inertíaalstelsel. Be- 


weegt S, eenparig versneld met versnelling a, t.o.v. S, dan is 


aj kk 
nek 
EL a, =ma -ma, =K-ma, 
â Willen we nu ook in zo'n stelsel S, de wet van Newton blijven toepassen, 


dan moeten we aan de kracht rd nog een fiíctieve- of schijnkracht (== ma, ) 


toevoegen. 








Ook in een stelsel, dat eenparig roteert t.o.v. een inertiaalstel- 


sel kunnen we de wet van Newton blijven toepassen, mits we schijnkrachten 


invoeren. 


Over de dynamica van één massapunt zullen we verder niet spreken. 
EL Wel nog in het kort íets over de dynamica van een stelsel vrije massa- 


punten. We onderscheiden de massapunten door een index Vv. Dan geldt: 


E | 


es . 
RE a 
Keke) Û Û Û 


dwz. we onderscheiden de kracht K, op het v-de massapunt in een uitwendig 


Ü deel Ke ‚ afkomstig van bv. uitwendige krachtvelden, en éen inwendig 


NEEN 





deel Ed afkomstig van de andere massapunten. 
En Td 
Dv vu 
He D 
waarin Ï,, de kracht ís door u uitgeoefend op v. 
Als Ke = 0 voor alle v spreken we van een afgesloten systeem. 
We nemen aan dat voor de inwendige krachten de wet aktie = reaktie zal 


gelden: 


ee = +8 


no Du 
en bovendien dat deze krachten langs de verbindingslijn v — u gericht 
zijn. 


We definieren vervolgens het zwaartepunt van het stelsel: 


_ 
mr pn 
dB mee Me) 
M v 


vo 


Uit de bewegingsvgl. volgt door sommatie: 


md 
vv v ‚Vv 
v v 


Û 
Maar 

) B = À a = O, zodat 

v v u | 
Veu 

> 
MR = j Et - ë 
Kb) 

Û 


Dit is de zg. zwaàrtepuntsstelling: Het zwaartepunt van een stelsel 
massapunten beweegt zo, alsof alle massa daar geconcentreerd was en alle 


uitwendige krachten daar aangrepen. (Merk op: als Ko afhangt van de 


aoe 


redde 
Zie € 


jn 


fn we 


ii 





coórdinaten E, moet men oppassen met deze stelling:). 


De totale impuls 
=> > 
PE a F = u. 
v vv 
v v 


Verder zien we dat: 


rd Pd u 
P = MR = rd (inpulswet) . 


Voor een afgesloten systeem ís derhalve: P = constant, dus R 5 at + b, 
waarin a en D willekeurige integratieconstanten. 

Kijken we ook nog even naar het impulsmoment. 
Het totale impulsmoment ís gedefinieerd als 


> > > 
L = 4 . 
) ov” Pp 


v 
Differentiatie geeft: 


. nd hd 
L= FAP. + FAD. = r 
Tp En ‚Te A 
v v 


Ke) 


de) 
tt 


>p u > i 
vo 


Û 


Laat zien dat op grond van onze onderstelling omtrent de inwendige 


krachten de laatste term nul is. Dus: 


Eu 

L = ) AS de (impulsmomentwet) . - 
‚Vv v 
KD) Kad 


Voor een afgesloten systeem ís derhalve het totale impulsmoment constant. 


Als we invoeren de coordinaten t.o.v. het zwaartepunt 


> > > > 
r.=r =R, zodat mr! = OQ 
v v Vv 


Û 


en dus 








dan definieren we het impulsmoment t.o.v. het zwaartepunt 


Bewijs nu zelf: 

Dan +ExP; 
in woorden: het totale impulsmoment om de oorsprong O is gelijk aan 
het totale impulsmoment om het zwaartepunt + het impulsmoment van de 
totale, in het zwaartepunt geconcentreerd gedachte massa, om O. 

Door naar t te differentiëren, vindt men (verifieer dit) 


Ei > 
L' ==) Far. 
Ee Kil 


v 
De impulsmomentstelling geldt dus ook t.o.v. het zwaartepunt, 
hoewel het coôrdinatenstelsel dat met het zwaartepunt meebeweegt geen 


inertiaalstelsel hoeft te zijn. 


Hoofdstuk I: Het princi van d'Alembert en de vergelijkingen van 


Lagrange. 


1. Het principe van d'Alembert. 


We bekijken weer een systeem van n massapunten, coordinaten B 
waarbij we nu toelaten dat de coordinaten niet onafhankelijk zijn. 
We nemen nl. aan dat er p ( p < 3n) betrekkingen zijn tussen de 


coordinaten, die verder ook nog wel van de tijd mogen afhangen: 


Fijgje ere KI En €) = 0. Jd = le..p (1) 


ET 


in OO 





„De, 


Voorbeelden: Een massapunt beweegt over een gegeven oppervlak (dat 
zelf ook nog wel mag bewegen), een slinger, 2 massapunten die een 
vaste onderlinge afstand hebben (halter). 

Men spreekt in dít geval van een holonoom systeem (rheonoom als 
de F's van t afhangen, skleronoom als dit niet het geval Ìs). Voor zulke 
systemen geldt dat het aantal vrijheidsgraden (d.i. het aantal coördinaten 
dat gegeven moet zijn om de configuratie vast te leggen) gelijk ís aan 
het wankel coordinaten dat onafhankelijk gevarieerd kan worden. In ons 
geval is dit Sn = pe | 

Bv. een bol op een volkomen glad plat vlak heeft 5 vrijheidsgraden 
(2 om het middelpunt vast te leggen, 2 hoeken om het raakpunt op de bol 
vast te leggen en | voor de rotatie om de verbindingslijn van middel- 
punt en raakpunt) en men kan deze 5 parameters onafhankelijk varieren. 
Het ís een holonoom systeem. 

Bij een bol op een volkomen ruw plat vlak, waarover hij alleen maar 
kan rollen is het aantal vrijheidsgraden nog 5, maar deze kunnen niet 
onafhankelijk gevarieerd worden. Onder de beperkende voorwaarden zijn 
er nl. niet-integreerbare betrekkingen tussen de differentialen van de 
parameters. Dit ís een voorbeeld van een níiet-holonoom systeem. Daar 
zullen we ons verder niet mee bezighouden. 

Fysisch gesproken worden de beperkende voorwaarden (constraints) 
gerealiseerd door reactiekrachten, die we in het algemeen niet kennen. 
Bv. de normale druk van een oppervlak, de spanning “in het koord van een 
slinger etc. Soms interesseert men zich alleen voor de beweging, soms 
ook voor de reactiekrachten. 

Bekijken we weet db wet vat Wentont 


—) 


ee A (2) 





== 


waarbij K, de bekende (in- of uitwendige) krachten voorstellen en 
KE de onbekende reactiekrachten. 

(1) en (2) geven ons 3n + p vergelijkingen, terwijl we ón 
onbekenden hebben Li en K) . Deze kunnen we niet oplossen. Het 
principe van d'Alembert zal blijken de ontbrekende informatie te 


leveren. 


Het luidt: Bij een kleine virtuele verplaatsing verrichten de 


reactiekrachten tezamen geen arbeid. 
Een virtuele verplaatsing Br, is een oneindig kleine verplaatsing 


op een vast tijdstip waarbij aan de voorwaarden (1) voldaan blijft. 

Er is geen algemeen geldend bewijs voor de juistheid van het 
principe van d'Alembert, men moet het opvatten als een onafhankelijk 
postulaat, dat door de resultaten gerechtvaardigd wordt. In een aantal 
concrete gevallen kan men de juistheid er van inzien: Hebben we bv. 
een bol op een glad plat wike De reactiekracht staat dan loodrecht 
op het vlak en verricht inderdaad geen arbeid bij een virtuele ver- 


plaatsing waarbij de bol op het vlak blijft. Verder: Stel we hebben 


2 massapunten op vaste onderlinge afstand (halter). De reactiekrachten. 


zijn hier langs de verbindingslijn gericht en gelijk en tegengesteld: 
ES pn ' 
K, en & id 
> 1 > > 
De door deze krachten verrichte virtuele arbeid is K, . (sr, - Br) 
—_} 


en dit is evenredig met {r; -r) ie B(r, -r) = is >, -r) = 0. 


Ook bij de slinger is direct in te zien dat de reactiekracht bij 


virtuele verplaatsing geen arbeid verricht. 


Het principe van d'Alembert luidt nu ín formule: 
) E'. sr, = 0 (5) 
P v v 
ze 


of 





=iÎ 


) (a, E, - E) e Br, = 0 (de) 


Ke) 


waarbij de er 's niet onafhankelijk zijn vanwege de beperkende voor= 
waarden. Merk op dat in (kb) geen reactiekrachten meer optreden. 
Sommerfeld geeft een aantal voorbeelden waarbij toepassing van (4) 
reeds leidt tot een oplossing van het probleem. 
3 
In de statica is r_ = O en dus ) Heer = 0. 
: v Vv v 
v 
Dit is het zg. principe van Bernoulli: In evenwicht verrichten de 


bekende gegeven krachten geen arbeid bíj een virtuele verplaatsing. 


Als de krachten afgeleid kunnen worden van een potentiaal 


dr, ae) 


KN Om en is) 
hb) hand 
8 


dan wordt het principe van Bernoulli: 


n 
Vv 
ar . BE =8V = 0 ‚ : (6) 
òr ie 
v=l ie 


Voor statisch*evenwicht heeft dus de totale potentiële energie een 
extreme waarde. 


2, De vergelijkingen van Lagrange van de 1° soort (1760) 


Waren alle Br, in (k) onafhankelijk dan zouden,we uit (k) kunnen 


besluiten tot 


rd 
nr = XK, 
vo Vv 


hetgeen niet juist is als er reactiekrachten zijn. We hebben echter 
nog nevencondities waaraan de virtuele verplaatsingen Br, moeten 


voldoen, ze zijn dus niet onafhankelijk. We maken nu gebruik van de 





wd 


multiplicatoren methode van Lagrange, 


Uit (1) r,- 0) volgt: 
af, òF òF. 
pn Re _l ot et Ì = 
oF, = }e grad F, Br, = i F, Bx, + A + 3, se) =O. (7) 
Û vv î 


Deze vergelijking vermenigvuldigen we met een voorlopig onbekende 


multiplicator \, en we tellen dan de ontstane betrekkingen op bij (hk). 


J 
n 
ng 2 zi 8 
} Re) n, Brad, F,) «ef, =O. (8) 
v=l j=1 
Er zijn Sn a p van de er 's onafhankelijk. We kiezen nu de Ass zodanig, 
dat p van de uitdrukkingen 
rd 2 5 
mrk,* De A, grad, F, = 0 (9) 


d=l 


worden. Uit (8) volgt dan dat ook de resterende 3n - p uitdrukkingen 
(9) nul moeten zijn omdat de bijbehorende (3n - p) Br 's wel onafhan- 


kelijk zijn. Zo hebben we 3n vgl. (9) gekregen die: de Lagrange vgl. van 


de 1° soort heten. Samen met F Fha O hebben we nu 3n + p betrekkingen 


tussen 3n + p onbekenden Ee, en X j 


We zien dat de reactiekrachten Juist gegeven worden door 
Ü ie f 
K, ) A, grad, F, {10) 
J : 


3n grootheden RK uitgedrukt in p A j 5 Het zullen in het algemeen 


functies van de tijd worden. 


Eek 


ii we 


Pe 


Eid 


ige 


En 


ke 


5 


ed 


En 


i bei 


ie eli 





nl 


RE be 


Toepassing: Stel er is een virtuele verplaatsing mogelijk (met in- 
achtneming van de bijvoorwaarden) waarbij alle punten dezelfde 


translatie ondergaan. Dus Br, = er. 


n 
eN 
Dus ) (e Rd ser = 0 voor iedere sr. 
vv Vv 
v=l 


Dan is 
n 
bard 
5 (a r -E) = © 
vv v 
v=l n 
py 
P = ) Ek 8 
v 
v=l 


Is de resultante van de bekende in- en uitwendige krachten nul, dan 


geldt dus Ë = constant. 


3. De vergelijkingen van Lagrange van de 2° soort. 


In de fysica zijn deze vergelijkingen, waaruit de reactiekrachten 
veraneni zijn, belangrijker. daureniaat men zich voor die reactie- 
krachten (spanningen) dan moet men dus teruggrijpen naar de Lagrange- 
vergelijkingen van de 1° soort. | 
We hebben weer 3n coórdinaten en p betrekkingen. We voeren nu ín 
Sn - p zg. gegeneraliseerde coördinaten a, (k = 1,...3n-p), waartussen 
geen betrekkingen meer bestaan en die de configuratie vastleggen. Dus 


evenveel gegeneraliseerde coördinaten als het systeem vrijheidsgraden 


heeft. Bij een vast lichaam bv. de 3 rechthoekige coördinaten van het 


zwaartepunt en 3 hoeken van Euler. 
â, heet gegeneraliseerde snelheid. De oorspronkelijke coördinaten 


Ky Yy Zo zijn functies van a, en eventueel van t. Dit laatste ís bv. 


het geval als de betrekkingen de tijd bevatten, maar ook als er geen 





- ik — 


betrekkingen zijn en als we bv. de gegeneraliseerde coórdinaten 


nemen t.o.v. een roterend assenkruis. De relaties FP, =O zijn 


J 
identiek vervuld ín de gegeneraliseerde coördinaten 





F, {a ae, yv, (apst)o eve (apt), «) = 0 
en dus 
n >p 
òr 
zi (eraa, F, 5?) = 0. 
v=l 


We vermenigvuldigen nu de vgl. van. Lagrange van de bel soort: 


Be -E, + a ); grad, F = 0 
j=Î 
= 


scalair met Se en tellen ze op (sommeren over v). 
0 


Op grond van (12) vallen dan de )\'s weg en we hebben 


> òr r 
\ GE E35) 

vv da v da, 
v=l 





voor 4 = Î,2,....3n =p = m. 
De 4 

Als nu FE (a, at) 

dan is m ze n 
=> àr ar 
Er) Ute 

v dq, 2 òt 
p=]l 


een lineaire functie van de â's- 


Hieruit volgt: 


> 
. 





(11) 


(12) 


(13) 


(ie) 


(15) 


























Dan is 
zi > ” je 
Een Dee (À =e) : ni) 
v àd v° dt \v° …e v° dt ° 
f) dâ, dd, dâ, 
Fy > 
òr dr 
ESC 2) En ze) (16) 
| > dà, ie 8 | 
Verder: 
> 25 m 2 
Ln) „) ha; 
dt da, da dt dat, k 
k=| 
e m > 
ar ar > 
EE. NE fak ee à _ Ò … 
da, \3t * da, à) 7 a, “o HE 
LA k=1 
(17) in (16) ingevuld geeft: 
> pg 
— àr > > 
en 
re N „re B 18) 
v a, dt \v di, v a, v 
(13) wordt nu: 
n sE n ar n sr 
d /e dd © ke k To 
) a, at Fo * ie Ken ) ken n (19) 
v=l £. v=i v=l 


n 
>è 
Nu is de kinetische energie T -) 2m n . Daardoor wordt (19) 
=Ì 


je | 


. e (20) 


a(aT)\_ AF È 

dt dq, % 
dq, ê 

Opmerking: omdat FE een lineaire functie is van de a's wordt T een 

kwadratische functie van de a's en wel als r (ast) niet expliciet 


van t afhangt een homogeen kwadratische. 








aib 


We schrijven nu het R.L. van (20) een beetje anders. De arbeid 


bij de virtuele verplaatsing Br, is 
n m 
jp 
BA= ) KR -8r, = ENC En) oa, +) aen 
v=l v=lg=l p=Î 


waarbij we 
n pd 
__Ù 
) Kk Er | (21) 
$ 
v=l 


de gegeneraliseerde kracht noemen. 


dr, 
Opmerking: Vergelijking van (1h) en (21) (als we afzien van <-” SE D) leert 


dat de transformaties van snelheden en krachten in de gegeneraliseerde 
snelheden en krachten samenhangen. Zij heten wel contragredient 
(= gespiegeld invers). 


(20) wordt nu: 


aft. B (22) 
de \ da, % 


de vergelijkingen van Lagrange van de 2°. soort. 

Uit deze Sn - p gewone differentiaalvergelijkingen van de 2° orde, 

die we afgeleid hebben voor holonome brite waarvoor d'Alembert geldt 

kunnen we in principe de a,(t) oplossen. Men kan (22) ook direct uit 

d'Alembert afleiden zonder via die van de 15 soort te gaan (vgl.Goldstein). 
Stel nu de krachten hebben een potentiaal, die nog wel van t mag 


afhangen 


É, ze grad, Vv Cajoeeezpt) « 


We kunnen de potentiaal ook in de a, s uitdrukken (vlajt). 


Dan is: 


a or BE 
EE 2 "9 AV 
A= À Kn s je ike he da, (23) 
v=l so 


ee 


ie 


Eni 





„iT 


Voeren we nu íinL=T =-Ve= Llaá,t) de functie van Lagrange. 


Dan kunnen we schrijven (daar V niet van â, afhangt ) 


ä „. %& | 
2E) ide (al) 
De oplossingen zijn bepaald door de beginwaarden van qa, en ä, 
De voorwaarden waaronder we de vorm (2) van de vel. van Lagrange van 
de 2° soort hebben afgeleid zijn: 
a) het systeem is holonoom 
b) het principe van d'Alembert geldt 
c) er is een potentiaal waarin de â, niet expliciet voorkomen. 


We definiëren nu de gegeneraliseerde impulsen: 


„Lb _ AT 
Pp, = = pn (25) 
k dd, A, 
Hiervoor geldt: 
n Fé n ag 
òr ar 
nE-) Bt =(volgens 15) ) EK. (26) 
y= Û v=l D 


We zien (vgl. 21) dat de krachten en de impulsen cogredient transformeren. 


Toepassingen: 


hi Een vrij deeltje in een conservatief krachtveld. Hier kan men bv. 


als gegeneraliseerde coördinaten de gewone rechthoekige coördinaten 
AL 


nemen. Dan is p‚ = in mk, de gewone impuls. 


2 e \ > À 
L = àm(% jes) = Víx‚y,z) 


De vgl. (2) worden hier: mx + N . 0 enz. dwz. de gewone vgl. van Newton. 


òx 





… 18 -— 


20 Neem in het platte vlak een vrij deeltje in een conservatief kracht- 


veld. Als gegeneraliseerde coördinaten nemen we vlakke poolcoördinaten 
r en p. 

2,2 2, 22 ' SRE ee 
L= Am(v‚ + ie Jm(t +rÒ) P, smf Po * 3 = mr} = mv, 


impulsmoment. 


De Lagrangevergelijkingen worden: 
mi - rj) + Lo 


òr 
m& (r°6) +S-o. 


Wat zijn de gegeneraliseerde krachten ín dit geval? 


NVN _PrELKre av 
Ce = da, R ë da, K da, + K, ë da, 


x= r cos @p y= r sin gp 


v 
-- Lex, 
ie 
"30" ko 


waarin K‚ en Ko de componenten ven de kracht in de r en de p richting 
voorstellen. o is dus eigenlijk het moment van de kracht om de oorsprong. 
Onze Lagrangevgl. zijn equivalent met de vroeger afgeleide bewegingsvgl. 
in poolcoordinaten. Als V = V(r) (centraal krachtveld) ís rò = constant, 


dwz. het impulsmoment mrd = constant. 


35° Leidt zelf de bewegingsvgl. af voor één massapunt in 3 dimensies in 


poolcoordinaten. 


h° Leidt zelf de bewegingsvgl. af voor 2 massapunten door een starre 
verbindingsstang zonder gewicht ter lengte # verbonden. Voer in de 
coordinaten van het zwaartepunt en de hoeken 6,p die de orientatie van 


de stang bepalen. 


KE Sea penn 


= (9 » 


à M (e tE + #2) + }I (e° + sin“e ) 





a _ 
m,m, 
12 @ 
M = m, + KN I= a, +, 8 - 


5? Roterend assenkruis 
We nemen de rechthoekige coördinaten 


xii in het eenparig roterende assen- 





kruis als gegeneraliseerde coördinaten . 


x xXx 


‚ cos wt =— y, sin wt De transformatieformules hangen hier 


x 


y sin wt + y, cos wt expliciet van de tijd af. 


1 





Te Am (is) Amir eum, - 2u te ty) 
T m T ° 
kaas sd en m(*, -uy) bi, 5 5, = my, tx) 


De Lagrangevergelijkingen luiden: 


. je B n AV 
mx, emuy, - Mw x, = òx, 
bed KJ è An € JV 
my, + emux, mwY, = àv, 


Hier komen de Corioliskracht en de centrifugale kracht te voorschijn 


praktisch zonder rekenwerk. 


hk. De energieintegraal en de hamiltonfunctie; cyclische coordinaten. 


Uit de. vgl. van Lagrange van de 2° soort volgt: 








Immers het linker lid luidt uitgeschreven: 


Ver Tak) -D Heder 
k k 


k „ & 


Op grond van de vgl. van Lagrange wordt dit - E. 


Als nu L niet expliciet van t afhangt (als bv. T noch V expliciet van 
t afhangen) dan vinden we: 


) &, ä, „LL = Je Ps â „L = Kla á) = constant . (28) 


k k 


Dit is de zg. bewegingsintegraal van Jacobi. 

In het belangrijke geval dat de transformatieformules van E naar 
de A 8 niet expliciet van t afhangen, wordt T homogeen kwadratisch in 
de de (vgl. opmerking p.16). Dan geldt volgens een stelling van Euler 

) dr. 
k 


In dit geval is 
) Pi â, „L=T+V (29) 


en is dus H gelijk aan de totale energie. 
Als bovendien e = O0 dan wordt de bewegingsintegraal van Jacobi juist 


de energieintegraal en hebben we 
Klad) = T + V = constant. (30) 
Het is gebruikelijk om door middel van Pk ke, ee anelieden in de 
functie Ì te elimineren en uit te drukken in de Py” Dan wordt 


Deze functie H noemt men de hamiltoniaan. Onder de genoemde voorwaarden 


a ET 


stelt H de totale energie voor van het systeem en blijft gedurende de 
beweging constant. (L niet expliciet van t af en T homogeen kwadratisch 


in de ä) 


Opmerking: Als een zekere coordinaat q,‚ niet in L voorkomt (ä, mag wel) 


Ê 
dan spreken we van een cyclische of ignorabele coordinaat. 
Dan ís U O en dus % =p, = constant. 
òQ, , 4, ' 
Voorbeeld: Stel a, is een cyclische coordinaat zodanig dat als a, met 
een zeker bedrag @ toeneemt alle xs met dat bedrag toenemen en de 
Y, 8e onveranderd blijven. Een verandering van a, stelt dus een 
translatie van het hele systeem ín de x-richting voor. De neven- 
condities moeten dit dus toelaten. Cyclisch betekent bovendien dat 
dan de potentiele energie niet verandert. De kinetische energie hangt 


nl. niet af van zodanige coördinaat. 


dE Fe Dr òx ày 
EL) (£ Ep Pr Se 
da, m, FE, ë aa, en da, = at Ja, = 0 want Ja, = Ì Ja, = 0 
KO 


Òz 
—_Ì = 0. Dus ar = o) e 
da, òd, 


a, is bv. de x coördinaat van het zwaartepunt, de andere q's zijn 


relatieve coordinaten. Nu is: 
> 
Ed or, sd dr, ° 
D, = nr az) = n(P SP) =) mk, = const. 
1 ‚_v\v” A4, v\v ag, ‚5 9 
* 8 “5 


Als dus een mechanisch systeem in zijn geheel een translatie in een 
bepaalde richting toelaat zonder dat de potentiële energie daarbij 

verandert, den is de component van de totale impuls in die richting 
constant. We zeggen dat het srate Anvartant is voor translatie ín 


die richting. Uit zo'n invariantie volgt dus een behoudswet. 








Be 


Iets dergelijks geldt als een systeem een rotatie om zekere as 
toelaat zonder dat de potentiële energie daarbij verandert. Invarian= 
tie voor rotatie. Dan is de component van het impulsmoment langs die 


richting constant. 


5, Elektron in een elektromagnetisch veld. 


De krachten hangen hier van de snelheid af; denk maar aan de 


Lorentzkracht. Als we echter een functie Vla ä” t) hebben zodanig dat 
a ES | (52) 
nt: at \ò, 
dan is het duidelijk dat uit 
LE 
at \d/ TR, ° 


volgt, als we weer invoeren L = T « V, 


ei) Ere 


d.w.z, de zelfde vorm van de Lagrangevgl. van de 2° soort, maar nu voor 


een gegeneraliseerd geval. We noemen 


ee _Àr _ àV 


nu de gegeneraliseerde impuls. Als L niet expliciet van de tijd afhangt 


is 

) Pi á, „-L =H 
weer een Ste H is echter niet meer gelijk aan T + V. Wel kan men 
laten zien dat H de totale energie voorstelt en niet T + V. (Kijk naar 
arbeid door de krachten Q k verricht.) 


Het hier beschouwde geval doet zich voor bij een elektron in een 


elektromagnetisch veld. 


PR: 





= 25 & 


We nemen voor 
> rde, 
VE, É, t) =e Jorr) - EAratd, (%) 
(e is de lading van het elektron), 


waarbij p en A de scalaire resp. vectorpotentiaal voorstellen. 


1 NR 
Neem maar: (Bedenk dat B = rot Â E - - grad P= ps He). } 
dA DA ÒA òA 
_.V, A/V) wek (er at 0 Mad 
A= VEL) - “xk edt em TE vt) 
òA DA dA dA DA DA SA 
DN Dek MET ee nde Eed. SOE WOE Dai EN a 
ek: ahd ad En mi 


u 
® 
Loe 
+ 
| 
<ì 
> 
dl 
Ld 
han, mnd 
hd 
, 
6. | 
ed 


en dit is juist de x-component van de Lorentzkracht. 
Ook voor een elektron in een elektromagnetisch veld gelden dus de 


Lagrange vgl. van de 2° soort met de lagrangiaan: 
TO 
Dehaes?) ef; - EE) (36) 


De gegeneraliseerde impuls 


jn mt e IT ( 
premie ed. (27) 
me >è jj 
De hamiltoniaan wordt: H=dmt +ep(r,t) 


è 
Ì 
ne {P- 2e, 0} E ep(F,t) „ TeV.” 





Hoofdstuk II Kleine trillingen. 


1. Bewegingsvergelijkingen en seculaire vergelijking. 


Eerst wil ik nog even herinneren aan ééndimensionale trillingspro- 
blemen die u reeds bent tegengekomen, bv. de harmonische oscillator of 
de mathematische slinger. Als we krachten hebben die alleen van x af- 


hangen dan hebben we steeds ook een potentiaal zodat 


K(x) en _ dv(x) 


dx . 


We spreken van een evenwichtsstand als er een waarde van x is waarvoor 


K(x) =O = ar ‚ We kiezen die evenwichtsstand als nulpunt en ontwikkelen: 
dv afv 2 
v(x) = v(o) + == . X +} 5 « + vec ce 
dx è 
X=O0 äx X=0 


Vo is willekeurig en kunnen we O kiezen. De lineaire term verdwijnt 
omdat we ontwikkelen rondom een evenwichtsstand. We beperken ons tot 


kleine uitwijkingen en breken af na de kwadratische term. 


Dus v(x) =Jax° (5) 
=Ó 


De bewegingsvgl. luidt dan: mx + ax = 0. 


Als V een minimum heeft in O, dan is a > O stabiel evenwicht. 


[a 
xx, sin ( Eiso). 


Als V een maximum heeft in O, is a < O labiel evenwicht. Exponentiele 


oplossing. 


We nemen nu een probleem met n vrijheidsgraden. De coordinaten hoeven 
geen rechthoekige cartesiaanse te zijn, we nemen gegeneraliseerde coordi- 
ten. Bv. dubbele slinger : hoeken Q en B. Verder toepassingen bij molecuul= 


trillingen, kristalroostertrillingen etc. We zullen zien dat we ook in 





=De 


het algemene geval van kleine trillingen het probleem 
kunnen terugbrengen tot dat van n ongekoppelde harmonische 
| p oscillatoren. | 
Onderstel: aq) Er is een potentiaal die niet van t afhangt, Vla,---a)- 
B) De functies Fat) hangen niet expliciet van t af. 
Dan wordt zoals we gezien hebben de kinetische energie 
homogeen kwadratisch in de á- 
T = Tav á)- 
y) Stel er ís een evenwichtsconfiguratie, waar de krachten 


nul zijn, dus 


2) = 0 voor alle k. 
En 


We geven de coordinaten En in deze configuratie de waarde nul. 
We gaan nu V weer rondom de evenwichtsstand ontwikkelen en breken af 


na de kwadratische term: 


vla,---a) =} ) bs A Ca) 
kt 
waar 
È 
Sy 
b = mm, = b pe (2) 


De kinetische energie is algemeen: 
Tek) og (Aert) de dg (3) 
k,@ : 
We kunnen ons hier de functies aeg djr An) ook weer in Taylorrseksen 
ontwikkeld denken en we zullen hierin alleen de eerste (constante) term 
meenemen, zodat we bij T evénver gaan (kwadratisch) als bij Vv. 


Ted) os ied hi 


kg 








waarbij we kiezen ag = ark . 
' 2e ze 

Als de q's Cartesiaanse coordinaten zouden zijn, zou ang ed mg en is 

(k) geen benadering voor kleine uitwijkingen maar exact. 


De Lagrangevergelijkingen van de 2° soort luiden: 
} ag a, + ) big de = Oe k = l,2...n (5) 
8 8 


Deze gekoppelde vergelijkingen kunnen worden opgelost door te stellen: 
1 
a,(t) =A, e“* (6) 
waarin de constanten w en A (op een gemeenschappelijke factor na) 
bepaald kunnen worden door invullen in (5). Onze coordinaten zijn 
uiteraard reëel. We nemen derhalve van de verkregen oplossing het reële 


deel (dit is zelf ook een oplossing). 


Na invulling vinden we: 
2 
NE), A, = O0. (7) 
p 


Hier staan n homogene lineaire vgl. met n onbekenden Ap die slechts 


een niet-triviale oplossing hebben als de determinant = 0 


act | - Pag tby| = 0 (8) 


de zg. seculaire vergelijking. Dit is een n-de graadsvergelijking voor 
a Noemen we a = \. We krijgen in het algemeen n wortels A, (J=l..en) 


en voor iedere j vinden we ook een stel A, . Als de,wortels verschillend 


Je 
zijn kunnen we de Asp bepalen (voor vaste j) op een constante factor ra. 


We noemen )\, de eigenwaarden en de rijen A 


; 30 de eigenvectoren van ons 


probleem. Voor iedere j hebben we een zg. eigentrilling of normaaltrilling, 


ook wel hoofdtrilling. 





„BT e 


2. Eigentrillingen en eigenwaarden. 


In matrixvorm geschreven luidt (7): 
anA = bA. (9) 


Als a de eenheidsmatrix was (as = Bg) dan stond hier het eigenwaarden- 
probleem van een symmetrische n bij n matrix. Ons probleem ís dus een 


beetje gecompliceerder.Vvoor een bepaalde j wordt (7) nu: 


) len Ns Dig) Azg =O (10) 
9 


(j) 
Noemen we ag A, - bg even Cp . Omdat de determinant van dit stelsel 
O is mogen we één van de vgl. wel weglaten. De Ass verhouden zich (j vast) 
als de minoren van de elementen van een willekeurige rij. Noemen we de 


minor van 6d) cd) dan wordt: 


ke k4 
A = C ctd) met 1 willekeurig (11) 
jk Jd ik i 
De en zijn reéel, want de ), zijn reeel zoals we direct zullen bewijzen. 
C, is een willekeurig complex getal. 
De j-de eigentrilling is nu: : 
iw,t 
EIN Aan (3) 
a = Re Anr e = |c| C‚x eos (w‚t + 9) 8 (12) 


9, is het argument van Ce ws, = SN: In een eigentrilling hebben alle 
a's dezelfde fase, dezelfde frequentie en een vaste amplitude verhouding. 
Ze gaan bv. tegelijk door de evenwichtsstand. De amplitude verhoudingen 
worden bepaald door de traagheden (eg) en de stijfheden (big) van het 
systeem. De algemene oplossing van ons trillingsprobleem is nu een wille= 


keurige superpositie van eigentrillingen. 


= 28 … 


Bewijs dat de A reéel zijn. Stelling van Sylvester: Als één der 


kwadratische vormen positief definiet ís, zijn de eigenwaarden X\ 3 reéel 


en zijn de eigenvectoren “orthogonaal". In ons geval is inderdaad T een 
positief definiete kwadratische vorm. 


Er geldt: 


* 
A; ) eg sg 7 ) beg Asg” vermenigvuldig met A,, en sommeer over k 
8 


8 


Ed had & 
Nn } as Ask = IR beg Air ‚ vermenigvuldig met Ap en sommeer over @ . 
k k 


Als we vervolgens aftrekken vinden we 
« « 
Oy A) ) ag Air Asg= 0 (13) 
kg 


Nemen we nu í = j en stellen we A Ask * iB, dan wordt de dubbele 


ik” AH 


som in (13): 


5 Eg % % * 8 Pr Pe tt } ag Pp TP) 
kg kg kg 


De laatste term ís nul vanwege de symmetrie van ag . 
De eerste beide sommen zijn nu positief definiet Z Eng ä, ä, > O voor 


â, #0. We zien dus dat uit (13) volgt: 


EN ie 


1 dwz. de eigenwaarden zijn reeel. 


We hebben zojuist gezien dat À as AL positief definiet is, 


ix San 
kf 


als jn ap Kk x, positief definiet is. 


* ad 
Uit A4 ) eg Aix Mag = ) big Aix Arg 
kg kf 

















volgt 
* ; 
VA Hip) ’ 
Een meat. (1) 
Tr Ar” Ag ) 
Is nu behalve de kinetische energie ook de potentiële energie nog 
positief definiet dan zijn alle A4 > O en dus alle w; reeel. In dit 
geval heeft de potentiële energie een minimum in de evenwichtsstand, we 
spreken van stabiel evenwicht en de oplossingen zijn werkelijk oscilla-= 
torisch. 


Uit (13) volgt verder voor Aj \; : 


ks 
he han hi = 0. (15) 
k4 


Dwz. de eigenvectoren behorende bij verschillende eigenwaarden zijn 
“orthogonaal" in de bovenstaande zin. We kunnen de Air nog zo normeren, 


dat we kunnen schrijven 


« 
ke ag Aix Asp = Dig (16) | 
kg | 

Uit As D ag As, = ) big As, volgt door vermenigvuldiging met 
0 4 
dr ‚ sommatie over k en toepassing van (16) 
DB Da Ar, Are (17) 
d Ad k4 ik 54 
k4 


De hier gegeven oplossing van het probleem der kleine trillingen is 
reeds afkomstig van Lagrange (1762). Heeft de seculaire vgl. (8) meer- 
en wortels, dan spreekt men van ontaarding. Lagrange dacht dat er 
dan oplossingen waren van de vorm Pelwt (zoals bij ééndimensionale 
gedempte slinger in grensgeval van aperiodiciteit). Pas in 1858 werd 
door Weierstrass de juiste oplossing gegeven: Stel dat \ 3 de enige dub- 


bele wortel is. Dan krijgt het stelsel lineaire vgl. de rang n-2 of 





m.a.w. we mogen nu twee A k willekeurig kiezen. We krijgen dan toch 


d 
het juiste aantal integratieconstanten. 


\, = O betekent dat deze trillingswijze geen eigenlijke trilling is. 


J 
Nultrilling, nulfrequentie. Men spreekt hier van indífferent evenwicht. 


Voorbeeld bij molecuultrillingen. De potentiële energie verandert niet 
bij braietente van het molecuul ín zijn geheel. Dit geeft aanleiding 
tot 3 nulfrequenties. Hetzelfde geldt voor rotaties. Bij een lineair 
molecuul geeft dit bovendien 2, bij een niet-lineair molecuul 3 nulfre- 


quenties. 


3, Hoof dassentransformatie. 


We kunnen het probleem der kleine trillingen ook aanpakken met behulp 
van een zg. hoofdassentransformatie, die op vele gebieden der natuurkunde 
(ook in de meetkunde) een belangrijke rol speelt. Men kan ook spreken 
van diagonalisatie van matrices. 

We zullen nl. laten zien dat we beide kwadratische vormen T en V 
(waarbij T positief definiet ís) door één lineaire transformatie gelijk- 
tijdig op hoofdassen kunnen brengen. Dwz. dat ín de nieuwe variabelen 


alle mengproducten verdwenen zijn. In de nieuwe variabelen worden dan de 


Lagrangevergelijkingen bijzonder eenvoudig, ze zijn nl.niet langer gekoppeld. 


Elke nieuwe variabele stelt een eigentrilling voor. 


Ik wil u eerst even herinneren aan het geval van één kwadratische vorm 


E bs Xe Xx, 2 


Men kan deze door een orthogonale transformatie 


Xe 5 J Pik 84 met 8 Pix Pik = 84; (18) 
1 k 


op hoofdassen brengen, dwz. dat we na uitvoering der transformatie 











'e 31 => 


krijgen : 

) At? | (19) 
i 21° 

i 


Bekijk daartoe het eigenwaarde probleem 


Ne, Ts NX, of ): (os - ADs) *4 = 0 en) 
£ > ú 8 


waartoe we komen als we vragen naar de extrema van je b 2 Xe X 2 onder bij- 


k/ 


B ehs (20) zijn weer n homogene lineaire vergelijkingen. 


voorwaarde je Xe 


k é 
Om ze op te lossen moet de determinant nul gesteld worden. We vinden n 


reële wortels i Bij elke eigenwaarde X i kunnen we dan een rij x-waarden 
vinden, die op een constante factor na bepaald zijn. We normeren deze 


eigenvectoren, die we p‚, noemen zodat we weer hebben (orthogonaliteit!) 
te | 


k 


We kiezen nu de matrix Pim die orthogonaal is, als transformatiematrix. 


Ei ) Pik OP KT ) Pix Br vel 
k A 


De eigenvectoren Pi Por’ .…… zijn de nieuwe basis. 


Nu hebben we voor de p's 


) bg P3g 5 Ma Pix 
P 


) bg Pik Pig * 3 ) Pik Pik * À\zbij- (aa) 
kf k 


Dan gaat de kwadratische vorm ) bg Xp over in: 





zE = 


T. 2 
TET, d Er 


kê i 


Opmerking: 
De lineaire vergelijkingen (20) hebben een eenvoudige meetkundige bete- 


kenis. De vgl. IN b, p XX Phs const. stelt een n-dimensionale ellípsoide 


kg 
voor. De uitdrukkingen Je bg x ri stellen de richting van de normaal in 


8 
het punt x, voor. We vragen in (20) dus naar die richtingen waarvoor de 


normaal juist langs de voerstraal valt. Dít ís iÍnderdaad de voorwaarde 
voor de hoofdassen van de ellipsoide. 

Nu hebben we in ons probleem te maken met 2 kwadratische vormen waar- 
van één positief definiet is. We willen deze door één lineaire (niet meer 
orthogonale) transformatie beide op hoofdassen brengen. 


Bekijk daartoe het eigenwaarde probleem 


NOME), xy = O0. (al) 
£ 


Merk op dat dit juist ons probleem (7) is. 
We noemen de eigenwaarden weer „, (reëel) en de eigenvectoren die daarbij 


horen p‚,_ ("orthogonaal'"). Normeer de eigenvectoren zo dat 
ik 


ä Seg Pax Pag T Pig» 


We hebben nu volgens (2) 


ke bg P3e * N ) eg P3o 
p P 


en hieruit volgt: 


D Dig Pik Pig * N ) Aeg Pax Pie * NM baje jaan) 
kg kg 








= SE 


Voeren we nu weer uit de transformatie 


* 5 ) Pix 5 (26) 
i 


dan wordt 


n 2 
ee 83 
kg 1 
b « A EE 
kg “Kk “0 ° ER 


ke i 


(27) 


Beide kwadratische vormen zijn op hoofdassen. 


Voor de nieuwe coordinaten E, worden de Lagrangevergelijkingen 


(a, - ) zet) 
Ym ki) 


i 


Si Hag > Da 


De &'s heten normaalcoördinaten, de assen ervan vallen langs de eigen- 


vectoren Pik Ik merk op dat wat we hier Pix genoemd hebben op pag. 28 


was Ask” 


A * ) Pik 4 7 Ee Pap Coon (uyt +9). 
i i 
Men kan de hoofdassentransformatie ook stapsgewijs uitvoeren. Eerst T 
door een orthogonale transformatie op hoofdassen. Alle eigenwaarden zijn 
positief. Dan gaan we een schaaltransformatie (contractie) toepassen en 
maken daarmee T tot een som van kwadraten. V is inmiddels in een andere 
kwadratische vorm overgegaan, die we nu door een niee orthogonale 


transformatie (die T invariant laat) op hoofdassen brengen. 





Eee ve 





Toepass ing: 


1e 


De dubbele slinger. 


We beschouwen kleine uitwijkingen uitslui- 


tend in het xy =vlak. 
2.2 

T, = m4 Ô, 
2 2 

T, = tm, (+95) À 


Nu is Xo = 6, sin 9, + 4, Sin Pp, 





Yo b, cos ®, + 8 cos Poe 


Y 
ER 22 p 
a è .@ 2 2 KJ . 
he” (4 % “4 Po + 24,40P, ô,) k 
2 2 
V= (m, *mo)g4,(1-cos p,) +m84,(1- cos p) =h(m, tm as, +Amo 4,0, 


Bepaal de eigenfrequenties. Laat zien dat N, en 5 positief zijn. 
Neem het speciale geval m, = m, b, = bo en bereken de eigenfrequenties 


en de amplitudeverhoudingen der beide eigentrillingen. 


\ 
Antwoord: pad = eef HS = 2 
® OO 
Ê AZ ee de JB, 

®, ®, 


m m m 


Gegeven: 5 gelijke puntmassa's (massa m), verbonden door 2 gelijke veren 
(veerconstante f). Beschouw alleen bewegingen in de richting van de 


verbindingslijn. Noem de uitwijkingen uit de evenwichtsstand respec- 
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tievelijk u u, ug Dan iss 
T= bm (ús + äs + ù5) 
vrktla sud. brlia) 
Beg Ss Bj gr Ups rn 


Bepaal de eigenfrequenties en de eigentrillingen. 


Hoofdstuk III Het principe van Hamilton en de Hamiltonvergelijkingen. 


1. Variatierekening. 


De tot nu toe behandelde mechanica was ín hoofdzaak reeds voor 1800 
bekend en gegroepeerd rondom de vergelijkingen van Newton en Lagrange. 
Omstreeks 1835 hebben Hamilton en Jacobi onafhankelijk van elkaar een 
algemene mathematische theorie opgesteld om de differentiaalvergelijkingen 
te integreren. Deze is als algemene theorie zeer belangrijk, doch wordt 
voor het oplossen van speciale eileen niet vaak toegepast. De theorie 
van Jacobi-Hamilton is ook voor latere onbielkkettng van de quantummechani-  *. 
ca van de arenie betekenis geweest. 

We hebben al eens gewerkt met een en differentiaalprincipe, n.l. dat 
van d'Alembert. Dit kwam ín de plaats van de grondvergelijking, terwijl 
de bewegingswetten er een gevolg van zijn. We behandelen nu enkele z.g. 
integraalprincipes, die direct de hele baan karakteriseren. We komen hier 
op het terrein van de variatierekening (Courant - Hilbert I, Frank = von 
Mises o.a. I, 881). Het grondprobleem is hier het volgende. We zoeken 


een functie y = y(x) ín een interval Ge,sXo) zodanig dat 
x 


2 
I= [ras ND) ax (1) 


Xi 


stationnair is. De functie F(x‚y,‚y') is bekend en yo), y(x) hebben 





voorgeschreven aen 

We leiden een noodzakelijke voorwaarde voor y(x) af. Zij y(x) de ge- 
zochte functie. Beschouwen we verder functies y(x) + an(x), waarin n(x) 
een willekeurige functie is, die nul ís in xj en x- De uitdrukking (1) 


wordt nu een functie van q. De noodzakelijke voorwaarde voor het extreem 


zijn van I is dan: 


GE), © 
je ok vof se 


Partiële integratie van de tweede term van de integrand; 


x 
ì {3-8 E,) nx) de = 0. 
1 


N 
Er geldt nu de volgende hulpstelling, die zich gemakkelijk laat bewij- 


zen. Als: B 


[ flx)nlx)dx = O voor alle mogelijke functies n(x), die vol- 


Kk 


doen aan núx,) = nx) = O0, dan zal f(x) = O0. Een noodzakelijk gevolg van 


het variatieprobleem 


5 


sf F(x,y,‚y') dx = 


Bi 


is dus dat y moet voldoen aan de zg. differentiaalvergelijking van Euler: 


dr AF _ o. 
dx ày' Ay 


Enige uitbreidingen : 


a, We kunnen vragen naar de van elkaar onafhankelijke functies 
a,(t),.-..-,a,(t), 


die de integraal 
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to 


I= dt F(t, ad) dt 


t, 


stationnair maken, waarbij a, (t‚) en a, (t‚) weer voorgeschreven zijn. 


Men schrijft Weers, 


2 
BI = s | F(t,aâ,) dt =O. 
t, 
Nu is: 
to 5 
… f- a : \ - 
I= ie ik {8 5a, + E. 5â, f dt = 0. 
t k 


1 


Partieel integreren geeft weer: 


t, | 
d 

ii ). EA - At B da,dt = 0. 

t, k 


Daar dit moet gelden voor willekeurige Ba, en de a onafhankelijk zijn, 


moet weer: 
s 
afd e E 
Se, iS, o DIE hs new (2) 


We zien hieruit, dat de Lagrangevergelijkingen feitelijk de Eulerverge- 


lijkingen zijn van het variatieprobleem 


sf rat = 0. 


Bij deze formulering van het bewijs zien we ook direct, dat omgekeerd, 


wanneer de differentiaalvergelijkingen van Euler (2) gelden, dan ook 
t ä 


2 
sf Fat = O0. 


8 


b. Hetzelfde probleeii met 2 onafhankelijke variabelen t'en x luidt: zoek 


de functie y(t,‚x) waarvoor 
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eff F(t‚x,Ys Yo Vig) dtdx = O0. (3) 
G 


De functie y(t,‚x) de op de rand van G voorgeschreven. In plaats 
van het discrete aantal vrijheidsgraden Aje geeft x nu continu 
oneindig veel vrijheidsgraden. Dit geval doet zich voor wanneer 
van het systeem van massapunten, 
zoals getekend in de figuur, over=- 
gegaan wordt tot het limietgeval, 
waarin co veel massapunten een on- 
eindig kleine onderlinge afstand 
hebben. 
Uit deze limietovergang kunnen we de bewegingsvergelijking van de trillende 
snaar opstellen. Een tweedimensionaal geval treedt op bij een trillend mem- 
braan, een driedimensionaal bij een trillende luchtmassa. 
Op de reeds eerder beschreven manier vinden“we de vergelijking van 


Euler voor dit vraagstuk: 


òr „AAF AF 0. 
dt av, * de ay, ° dY 


xXx 


(b) 


De Lagrangiaan voor de trillende snaar ís: 


vf{se (B) «0 (E) } ee 


waarin het eerste deel de kinetische energie is en het tweede deel de 
potentiéle energie. (5 = massaf/ cm, S = spanning). Stellen we nu 


[rax = L, dan gaat (k) over in de golfvergelijking 


Az 


è è 
òt òx 


(5), die hiermee equivalent ís, wordt 


s frat = 0. 





Het probleem van de trillende snaar hoort feitelijk niet meer thuis in 
de klassieke mechanica — waarín slechts gewone differentiaalvergelijkingen 


optreden — maar ín de veldentheorie (partiële differentiaalvergelijkingen). 


c. Evenzo kan de theorie uitgebreid worden, wanneer F nog van hogere af= 


geleiden van de gezochte functie afhangt. (Frank - von Mises, p. 886). 


d. Ook kan het nog voorkomen dat de gezochte functies aan nevenvoorwaarden 
voldoen. Als voorbeeld kunnen we vragen naar functies x(t), y(t) en z(t) 


die voldoen aan 


to 


sf F(t‚x,y,z,k,9, 2) = 0 
t, 


met als nevenvoorwaarde: G(x,y,z,t) = O0. 


Dit laatste stelt een ín de ruimte bewegend oppervlak voor als x, y‚, en z 
de coordinaten zijn met de tijd t als parameter. Dit probleem lossen we 


op met de multiplicatoren-methode van Lagrange. We vinden weer: 
2 . 
[ &E-Ele{AB-Bjor{âE-Eleere 0 


Bx, Sy en 5z zijn voor vaste t niet meer willekeurig te kiezen. Ze moeten 


nl. voldoen aan: 


U 


B er + Lays Boz 0. 


Xx ày dz 
Vermenigvuldiging hiervan met \, gevolgd door integratie over t en op- 


telling bij (5) geeft: 
to 
a FP dE 5 
hi BES EnEl Benet dt = Os 


We kiezen \ zodanig dat de eerste term in de integrand gelijk wordt 


aan nul. Tengevolge van de onafhankelijkheid van &y en 8z, verdwijnen 











“z ho - 


dan ook de beide andere termen: 


d n EG _ 

HE EN 4 

d 

ri Ae ee (6) 
da … f GE _ 

zE rh 0: 


Opmerking: Generalisatie voor meer functies en meer nevenvoorwaarden 


ligt voor de hand. 


2. Het principe van Hamilton. 


Keren we nu terug tot de mechanica. We bekijken nog eens de verge= 
lijkingen van Lagrange van de eerste en díe van de tweede soort. 
Voor n massapunten met p relaties tussen de coördinaten hadden we: 
pn + z = . 
mr a | ), grad, F, o Flott) o 
j=l 
Voor een systeem met een potentiaal: K, ze grad, V is deze vergelijking 


ook te schrijven als: 


p 

d … L je 

zE ed zh ) A srad, F, = 0. 
d=l 


evenzo voor Z . 
f Sy ”, 


Deze uitdrukking hadden we ook gevonden als we de Euler vergelijkingen 
hadden opgespoord van het variatieprobleem met nevenvoorwaarden: 
t, E 
sf L(t,P SE) dt = O0. FI Zot) =0. 


t, 


(vergelijk (5)). 
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Hiermee is dus aangetoond dat de vergelijkingen van Lagrange van de 
banen soort niets anders zijn’ dan de Eulervergelijkingen van bovenstaand 
variatieprobleem. 
Voor een systeem met een L=functie is, als we de gegeneraliseerde coordi= 


naten A noemen, waartussen geen relaties meer bestaan: 


a % _ U sG 
dt, A, 
en dit zijn juist de Eulervergelijkingen (2) behorend bij: 
to 
8 be Ldt = O. (geen nevenvoorwaarden) 


t, 


We hebben hier het zg. principe van Hamilton (183) gevonden. Voor een 
holonoom systeem met een functie van Lagrange, volgen de vergelijkingen 
van Lagrange (voor al of niet onafhankelijke coórdinaten) dus direct uit 
het principe van Hamilton, als dit tenminste voor juist wordt gehouden. 


Ook omgekeerd volgt uit de vergelijkingen van Lagrange dit principe: 
t t 


t t 
2 2 2 2 
a n Ld } 
fra - á SL dt = ) Een + RE Ee Ba, dt. 
‘ k k k 
t, t, k be Hp S 





Na toepassing van de vergelijkingen van Lagrange blijkt nu, dat 
t 


md Pi Öl e (7) 


in) k he 


Verstaan we onder $ een zodanige variatie, dat daarbij de grenzen niet 


gevarieerd worden (sa, (t‚) = a, (t) = 0) dan vinden we dus: 
t ä 
8 [ L dt = O0. 


t, 
Door dit principe van Hamilton wordt de hele beweging gekarakteriseerd. 
Een voordeel is, dat in deze formulering geen bepaalde keuze van coór= 


dinaten is gedaan. 





… h2 … 


tt 
2 
De uitdrukking hi L dt kan opgevat worden als een functie van begin= 


t, 


coördinaten, eindcoördinaten en van het tijdsverschil: 
to 


1, 


In deze variabelen uitgedrukt heet de functie W de principal function 
van Hamilton. 

Er zijn nog een abel andere integraalprincipes, die alle gemeen 
hebben dat de te variëren grootheid de dimensie energie maal tijd = actie 
heeft. (In't Duits heet dit "Wirkung'"; ook de constante van Planck h 
heeft deze dimensie "Wirkungsquantum'".) Men vat deze integraalprincipes 
wel samen onder de naam principes van de kleinste actie. Dat we hier wer- 
kelijk met een minimum te maken zouden hebben, is niet altijd het geval. 
Een beschouwing hierover vindt men o.a. in Whittaker. Wij zullen liever 


è 
zeggen: 5 f L dt = O betekent dat de integraal stationnair is. 


t, 


We kunnen het principe van Haniiton enigszins aanschouwelijk maken 
door de n onafhankelijke coördinaten a, op te vatten als rechthoekige 
coördinaten van een n = dimensionale ruimte, de zg. configuratieruimte. 
De configuratie van het systeem op een bepaald tijdstip wordt gekarakte- 
riseerd door een punt ín deze ruimte. De beweging van het systeem wordt 
dan voorgesteld door een "baan", díe in deze ruimte door het representa= 
tieve punt wordt doorlopen. Laat het systeem nu eens bewegen van P‚{t‚) 
naar P‚(t‚). Deze beide punten leggen 2n integratieconstanten vast, dwz. 
ze karakteriseren de bea eenduidig. Naast de werkelijk doorlopen baan 
kunnen we gevarieerde banen beschouwen, die tussen akrmidde eindpunten 
en eveneens in hetzelfde tijdsinterval verlopen. Het principe van Hamilton 
zegt nu, dat als we de werkelijke baan met deze gevarieerde banen verge- 


2 
lijken, hi L dt voor de werkelijke baan stationnair is. 
t 
1 











zùs = 


We zullen nu het principe van Hamilton nog eens direct uit dat van 
d'Alembert afleiden. Dit laatste luidde: 
_— 
…. pe p je 8 
) (mr, Ek) Br, = 0, (8) 
v 
> 
waarin K, de bekende in- en uitwendige krachten zijn. Br, is hierin het 
verschil in Be voor werkelijke en gevarieerde baan op hetzelfde tijdstip 


en moet nog aan eventuele beperkende voorwaarden voldoen. Door de sub- 


stitutie 
> > > ee 
nr dd et) helt? 
Br, at (Et sr) Fo at Br, = at Sr) 58(E,) 


gaat (8) Oer in 


bn À >} 
d ° > EN „è pd nú 
At Pe = + e = + e 
AT ij mn ( g sr) jj sm 8, Je K_-br, = ST + 8A 
v 5 v v 
ST is het verschil in kinetische energie bij het doorlopen van de werke- 
lijke en van de gevarieerde baan op hetzelfde tijdstip. BA is de door de 
krachten K‚ verrichte arbeid bij de virtuele verplaatsing Br, . Integra- 


tie van t‚ > t, geeft 


1 2 
to dea t 

f (ST + 8A) dt = ) Boer, | % 
t, v t, 


Pd 


r > sj 
Omdat srt) = sr (to) = 0 wordt dít 


to 
je (ST + 8A) dt =O. 


t, | 
f er é 
Schrijven we nog SA = 4 Ì, . ör,, = en ed, ‚ waarin Ee de reeds vroeger 


gedefinieerde gegeneraliseerde kracht voorstelt, dan is 





_ bh 


to 


f cer+) ep) a -o0 

t, kk 

de meest algemene vorm van het principe van Hamilton, ook geldig als 
de krachten geen potentiaal bezitten. De Eulervergelijkingen voor dit 


probleem zijn 
srt 
S= OE Os 
ät Ab, a, % 
dus de algemene vergelijkingen van Lagrange van de 2° soort. Voor krach= 
ten die wel een potentiaal hebben ís SA = =- 8V, zodat het principe van 


Hamilton weer de bekende vorm krijgt: 


t, t, 
1 sS(T = V) dt =O of sf L dt = 0 


t, t, 


(t, en t‚ invariant). 


2 


We zijn nu uit het principe van d'Alembert via de Lagrangevergelijkingen 
tot het principe van Hamilton gekomen, en ook uit d'Alembert via Hamilton 


tot de vergelijkingen van Lagrange. 


3, De vergelijkingen van Hamilton 


We gaan uit van de vergelijkingen van Lagrange 2 soort, n 2° orde, 
gewone differentiaalvergelijkingen. Deze kunnen we meer algemeen schrijven 


als 
f‚(t‚aâd) = 0, 1 = Is eeerN e 


Een dergelijk stelsel vergelijkingen kan altijd vervangen worden door 


een stelsel van 2n vergelijkingen van de 1 orde. Als we nl. ä, = Yx 


als nieuwe afhankelijke variabelen invoeren, krijgen we 





_h5 = 


U Ie 


In de mechanica kan dit nu bijzonder mooi gebeuren, door naast de a 
niet de ä, maar de P, als nieuwe variabelen in te voeren. 


We hebben L = L(q,,4,,t) en we voeren ín p, = er ‚ Uit deze vergelijkingen 
k’ “k k 


od, 
denken we ons â, opgelost, zodat we dan 


k 
kunnen beschouwen als een functie van A Py Cn t. Deze functie H(p At) 


heet de Hamiltoniaan of de Hamilton-functie. Dat de á, uit FP * % kan 
k 

worden opgelost, kunnen we gemakkelijk laten zien ín het eenvoudige geval, 
dat V niet afhangt van de snelheden en T homogeen kwadratisch is in de 
snelheden. Immers dan is Pk * 5 en de ä zijn oplosbaar, als de deter= 

er k 
minant sar! „ O0. Dit is nu inderdaad het geval, want T is een posi- 

k "6 
tief definiete kwadratische vorm in de snelheden en een noodzakelijk ge- 
volg daarvan ís dat genoemde determinant positief is. 
In het geval dat V niet van de snelheden afhangt en T homogeen kwadratisch 
is in de snelheden is de Hamilton functie H gelijk aan de totale energie. 
Immers: 

H = pâ, == L= ST 4 bete teValtY. 
‚Ek'k dd, k 
k k 


Altijd geldt: 


k k k k 
Nv jk ko 
= ) aap, ) da, da, Dt dt . 

k k 


Hieruit volgt nus 








à-À EX El, 
oP. da, ° a tt 
We zien hieruit dat als L niet expliciet van t afhangt dit ook voor H 
geldt en omgekeerd; verder, als de niet expliciet in L voorkomt (dus 
een cyclische coördinaat is), hangt ook H niet van A af en omgekeerd. 


Met behulp van de vergelijkingen van Lagrange die we kunnen schrijven: 


P = L ; 
k 
vinden we nu de zg. vergelijkingen van Hamilton, of kanonieke bewegings- 
vergelijkingen: 


5 be 8) 


Dit zijn nu 2n differentiaalvergelijkingen van de 1° orde met 2n onbekende 
functies. Ze hebben een bijzonder eenvoudige vorm omdat de differentiaal 
quotiënten al meteen zijn opgelost. 

Px En 4 noemt men wel kanoniek aan elkaar tousebesië variabelen. 

We zijn hier van de Legrangiaan op de Hamiltoniaan overgegaan door middel 
van een zg. Legendre-transformatie, die veel ín de natuurkunde optreedt. 
In zijn eenvoudigste vorm krijgen we zo'n transformatie als volgt: 


We hebben een functie f(x,y) van 2 variabelen x en y. We noemen nu: 


„af pn 
2 = == YZ Êf 


Uit z = & kunnen we y oplossen als functie van x en z, zodat we g kunnen 


opvatten als functie van x en z. 


dg = ydz + zdy - & dx = n dy = ydz - 8 dx 


zodat 


Oe OO 





u 


mimi d 


C 





-h7 - 
6 
ook de overgang van g(x,z) naar f(x,y) is weer een Legendre-transformatie. 
Het is duidelijk dat de overgang van Llaât) naar H(aoP»t) er een 
voorbeeld van ís. Bekende toepassingen vinden we in de thermodynamica: 
Van U(S,V) naar F(T,V) = U = TS, vandaar naar G(p,T) = F + pV. 
Behoud van energie. Bekijken we eens de afgeleide van de Hamilton- 


functie naar de tíijd: 
dH _ aH , aH H 
EE En) E. 
k k 
k : 
Als we hierin de Hamiltonvergelijkingen (10) invullen dan vinden we: 


aH _ A 
dt ° at ° 


Als H (of wat op het zelfde neerkomt L) niet expliciet van de tijd af= 
hangt, verandert H hiet gedurende de beweging, m.a.w. is H een integraal 
van de beweging. In het speciale gevel (maar veel voorkomend) dat 

H = T + V gaat dit over in de wet van behoud van energie. 

We hebben vroeger ingevoerd de n-dimensionale configuratieruimte. Een 
punt hierin geeft de configuratie van het systeem aan, maar zegt niets 
over de snelheden. Als de Lagrange=functie gegeven is, zijn door elk punt 
van de configuratieruimte nog oneindig veel banen mogelijk. Als we verder 
een baan in deze ruimte geven, weten we welke configuraties worden door- 
lopen, maar we weten nog niets over het tijdelijk verloop van de beweging. 
In het formalisme van Hamilton ligt het voor de hand om een 2n-dimensionale 
ruimte van de a's en de P's in te voeren, de zg. phaseruimte. Een punt 
hierin geeft zowel coördinaten als snelheden atie Hal phasgebaan karakteri- 
seert de beweging volkomen, geeft behalve de opeenvolgende configuraties 
ook voor elke configuratie de impulsen (en dus ook de snelheden) aan. Door 


elk punt van de phaseruimte gaat î en slechts | baan, immers 
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dg, = a dt dp, == el dt , 

k òP, k da, 
dus uitgaande van een gegeven punt geven deze vergelijkingen de raaklijn 
in dat punt. 

We kunnen eens de vraag opwerpen, hoe een willekeurige mechanische 


grootheid f (po dot) gedurende de beweging verandert. We vinden: 


eG En) 5 dt ” ,& En EEDE 


(11) 


Voor twee willekeurige functies van coordinaten en impulsen is de haakjes- 


uitdrukking van Poisson gedefinieerd als: 


(£‚g) = ) € 5 - zl E). | (12) 


Met behulp van deze schrijfwijze gaat (11) over in: 


ate, (15) 


en at 


(13) is de meest algemene vorm van de bewegingsvergelijkingen en leert 
Bis hoe een willekeurige mechanische grootheid met de tijd verandert. 
Als we voor f a, en P‚ nemen, krijgen we de Hamiltonvergelijkingen terug, 
ook f =H geeft een bekend resultaat. 

Het analogon in de quantummechanica van de Poissonhaakjes is de commu- 


tator: 


f = f hd f . 
Ì op?’ 6 | op Sop ” Fop “op 
Deze spelen in de'quantummechanica een fundamentele rol, ‘allerlei stel- 


lingen die we hier met Poissonhaakjes formuleren (bv. (13)) blijven in 


de quantummechanica geldig als we (f‚g) vervangen door  [f,e)- 





_ Voorwaarde voor een bewegingsintegraal is bv. commuteerbaarheid met 


de H=operator. 
Uit (13) volgt dat de noodzakelijke en voldoende voorwaarde opdat een 
functie f(po a) die niet expliciet van de tijd afhangt, een bewegings- 


integraal ís, dwz. constant blijft gedurende de beweging, luidt 
hard 
(É‚H) = O. 
Eigenschappen van Poissonhaakjes: 


Dit volgt direct uit de definitie (12). De quantummechanische analogieen 


zijn de fundamentele verwisselingsrelaties van Heisenberg. 


(f‚g) == (gf) 
(Ff + g‚h) = (f‚h) + (g,h) 
(fg‚h) = (f‚h)g + f(g, h) . 


Verder de zgn. identiteit van Jacobi (voor bewijs zie men Goldstein): 
(£‚ (g,h)) + (g, (h‚£)) + (h‚ (£‚,g)) = 0. 


Hieruit volgt onmiddellijk de stelling van Poisson: Als f en g beide 
bewegingsintegralen zijn, geldt dit ook voor (f‚g). 


Immers: 
d: 
at (f‚g) == ((£,g), H) hack ((g,H) ‚£) hk ((H‚£), g) = O0. 
Men late zien dat voor de componenten van het impulsmoment geldt: 


L, = (Lo Li), enz. In een centraal krachtveld zijn alle drie componenten 


bewegingsintegralen. 


Hoofdstuk IV De viriaalstelling. _ 


(vgl. ter Haar, Elements of Thermostatistics, n 


pag. 10. 
Becker, Theorie der Wärme, pag. 85). 


“Tenslotte zullen we nog kort bespreken de zg. viriaalstelling. Deze 
behoort eigenlijk al niet meer tot de mechanica, omdat het hier gaat om - 
gemiddelden. We bekijken een systeem van massapunten. ee 


De bewegingsvergelijkingen zijn: 
Eek rt Rt 
= = + . an 
Bo Fo v v v (1) 


We hebben de kracht gesplitst in een uitwendige kracht K\ (afkomstig van 
eventuele uitwendige krachtvelden of bv. ook van de wand van een vat) en 
een inwendige kracht É afkomstig van de onderlinge wisselwerking van En 


de massapunten. 


Nu is: Ee 
4 > > 1 > u == 
Je nr er = ) rr .K = ) r KS + F.R. (2) 
vv v v v v v NE ' v 
v v v v ed 
Verder is; B Oe 
> e —ê 
> . d Dad nd . ha 
m r . FT Sn mr , Y ee, mr ' 
je vv v dt ) vv v ) vv 
v v v Ì 
ê 
=$ 3 ) ar? „2T, (3) 7 
ät vv « a 
ö 
0 : —: 
waarin T = } ) mf, de kinetische energie voorstelt. ne, 


ov: f en 
We substitueren (3) ín (2) en gaan het resultaat middelen over een lang 


tijdsinterval. Het gemiddelde wordt aangeduid door < >. 











de 


t 


ODE n K, > heet het viriaal (ingevoerd door Clausius) 








KD) « 
t, 
d —2 1 ad 2 
en È ) mr,>= Lim EPS? fs 3) mr = 
dt kkn RA de ke 
t, 
1 d 2 
Se ee tot, (&3 22) | 
zie Wia v t, 





1 
ES Lim ir ( 


toet, > 00 2e 1 


In het algemeen zal de: uitdrukking ä mr 
nn 
eindige waarde hebben, als tenminste geen deeltjes naar het oneindige 


. on op elk tijdstip een 


gaan en de snelheden ook eindig blijven. 


In dat geval gaat de uitdrukking (kb) in de limiet naar O en vinden 
we: 


)r Bree) E ER d ee KK me del>. 
v v v v v v 


Û Û ke 


Dit heet de viriaalstelling. 
We zullen twee toepassingen bespreken: 
Dn Stel er is een potentiële energie V, die een homogene functie ís 


ad 


van de graad g, dwz. 


> > > n g 
var, Et), =z a vr, Roed) . 


Differentieer dít naar a en stel daarna a =| : 


t ( 8, 








- 52 = 


Derhelve: n B 
Viriaal = <) rk > e= =g<V> 
ke 


en de viriaalstelling gaat over in 


Kijken we bv. naar de Keplerbeweging in een ellips EF blijft nier eindig). 
De potentiaal is hier VEEN ee S, dus g ==. 
Dan geldt: 2 <T> + <V> = 0. 


Inderdaad volgt uit een gedetailleerde berekening: 


Pon 
p a 

…- @ 

<i> e= & 


en de (constante) totale energie is - EN . a ís hier de halve grote as 


van de ellips. 


2. Bekijk een gas (of vloeistof) bestaande uit N molekulen in een vat 
met volume V. De uitwendige krachten zijn alleen afkomstig van de wand 


en bedragen -p (druk) per oppervlakte-eenheid. Dan is 


áp E ° Re > =ep [ras (oppervlakte integraal) . 


Door toepassing van de stelling van Gauss volgt hieruit voor het uit- 


wendige viriaal: 


De viriaalstelling geeft dus: 


3pV=2<T> + <) ER: (7) 


v 
Nu is, zoals bekend (en we later nog zullen zien) 
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id ie. en 
Ae EN 
Y 


mn 
/ 


2e<T> = E NKT 


waarin T ín het rechter lid nu de absolute temperatuur voorstelt. 
Voor een ideaal gas zijn er geen inwendige krachten en daarvoor vinden 


we de ideale gaswet (wet van Boyle, Gay-Lussac) 


tm 


PV = NKT . (8) 


‘ 


Voor een werkelijk gas zijn er wel intermolekulaire krachten en kunnen 


ES 


we in beginsel met behulp van het inwendige viriaal correcties op de 


ideale gaswet berekenen. Men schrijft vaak: (v = k ‚ het volume per deeltje) 


ia En 


As 


Pv _ « . BT) , C(T) 
ek de de dake (9) 


B(T), C(T) .… heten de 2%°, zde 


er dn 


‚ etc. viriaalcoefficient . 


In Becker $29 worden de correcties op de ideale gaswet bij benadering 


berekend en wordt zo de toestandsvergelijking van Van der Waals afgeleid: 


(» + el (+-») = KT (10) 


waarin a en b voor het gas karakteristieke constanten zijn. 


We kunnen (10) ook schrijven: 








kT a 
DP y-b 2 
v 
of: | 
re aid 
p be 2 3 À 
j 
ij EE: kr È P+ 2 beden (11) | 
vo ov | 
À 


zodat voor de vergelijking van Van der Waals de viriaalcoefficienten 


zijns « 


| C=b, D=b", etc. 











| 
[ 
d 
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Werkkollege Statistische Mechanika 1974-1975 


(Vrije halter). Beschouw twee gelijke massapunten, door een starre verbindings- 
stang zonder gewicht ter lengte 1 verbonden, die zich verder vrij in de ruimte 
kinnen bewegen. 
Voer ín een geschikt stel gegeneral iseerde koördinaten (evenveel koördinaten 
als het systeem vrijheidsgraden heeft). 
Schrijf op voor dit systeem: 

De Lagrangiaân. 

De Lagrangevergelijkingen van de le en van de 2e soort. 

De Hamiltoniaan. 

De Hamiltonvergel1i jkingen. 

Welke koördinaten zijn cyklisch? 

Welke behoudswetten volgen hieruit? 


Vindt u zo alle behouden grootheden? (Aanwijzing: In een systeem met n 
vrijheidsgraden zijn er 2n-1 onafhankelijke behouden grootheden.) 


Bereken de kracht die de halter uitoefent op de twee massapunten als funktie 
van de tijd en de beginvoorwaarden. (Los daartoe op de vergelijking van 
Lagrange van de le soort.) 


Een dubbele slinger bestaat uit een mathemâtische slinger ‘lengte l,, massa m) 


1 , 

waaraan bevestigd is een tweede (lengte lo, massa m). Het ophangpunt van deze 

slinger is het massapunt van de eerste. 
Beide slingers bewegen in één plat vlak. 


Bepaal de eigenfrequenties van dit systeem. 


# 


Bereken voor kleine me de eigentrillíingswijzer van dit systeem. 
Bestaat er een gegeneraliseerd koördinatenstelsel, zodanig dat één van de 
twee koördinaten cyclisch ís? (Eventueel in bepaalde limietgevaller.) 


Aanwijzing: Bewijs dat alle Ed O. {Algemeen geldt dat een cyklische 


Ĳ 


koördinaat aanleiding geeft tot een eigenfrequentie w gelijk aan nul, en 


omgekeerd. ) 


(Roterend assenstelsel). 

Beschouw een vrij deeltje in een conservatief krachtenveld. Beschrijf de bewe- 
ging van het. deeltje (volgens Lagrange en Hamilton) ín een met eenparige snel- 
heid w draaiend rechthoekig koördinatenstelsel. 

Bereken zo de Corioliskracht ‘en de middelpuntvliedende kracht. 


Hoe kan men experimenteel een onderscheid maken tussen deze twee krachten? 
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Werkkollege Statistische Mechanica 1974-1975 (vervolg 3.) 


Veronderstel dat een waarnemer zich bevindt in een afgesloten ruimte 
(een "Skylab" zonder vanen bijvoorbeeld), die met eenparige snelheid 
w roteert om een vaste as, 

Bedenk een eenvoudig proef je hoe deze waarnemer zovel w als de plaats 
van de draaias kan bepalen. (Verschillende mogelijkheden. ) 


(Halter op 2 rails). 
Twee gelijke massapunten, verbonden door een halter van lengte 1 
bewegen zich zonder wrijving elk op een rail. De twee rails kruisen 
elkaar onder een rechte hoek op een afstand a (a < 1). De kortste 
verbindingslijn van de twee raíls is loodrecht op de richting van 
het zwaartekrachtsveld, terwijl de rails er elk een hoek van 45° mee 
maken. 
Beschrijf de beweging met de Lagrangevergelijkingen van de le en 
van de 2e soort, en met de Hamiltonvergeli jkingen. 
Bereken de kracht die de halter uitoefent op de twee massapunten, 
uitgedrukt in de gekozen gegeneraliseerde koördinaat. 
Met welk reeds eerder beschouwd systeem is dit systeem mathematisch 
equivalent? 
Bediscussieer de mogelijke bewegingswijzen van het systeem als funktie 
van de energie. 


Beschouw hetzelfde systeem als in opgave 4, echter nu met de verbindings 
stang van de halter vervangen door een veer met rustlengte 1 en veer- 
konstante w. Neem a = 0. 

Schrijf de Lagrangevergelijkingen van de tweede soort op. 














ze : Ed 
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Werkkollege Statistische Mechanica 1974-1975 (vervolg 5.) 


en los ze op | eed 
a) voor 1 = 0 ' 7 
b) als de versnelling van de zwaartekracht nul is. 


Bepaal voor kleine uitwijkingen uit een evenwichtsstand 





de eigenfrequentie (s) van het systeem. 





4 Drie gelijke puntmassa's zijn door vier gelijke veren 


verbonden met 2 wanden. Zij bevinden zich op één lijn, 








\ die loódrecht op die vanden staat, en 
kunnen alleen langs die lijn bewegen. 


Bereken de eigenfrequenties en eigen- 





trillingswijzen van dit systeem, en laat 
zien dat de laatste orthegonaal zijn. 
Onder welke omstandigheden is deze oplossing exact, ook 


voor grote uitwijkingen? 


7. In een homogeen zwaartekrachtveld bevinden zich twee 





punten op niet gelijke hoogte, die zijn verbonden door. 
een gootje, waarin een puntmassa wrijvingsloos kan glij- 
den. De vorm van het gootje is dusdanig dat de valtijd 
van de puntmassa vanaf het hoogste punt naar het laagste 
minimaal is als de puntmassa zonder snelheid wordt los- 
gelaten in het hoogste punt. De curve die het gootje be- 
schrijft heeft een brachistochroon. 

Stel de differentiaalvergelijking waaraan de 

brachistochroon voldoet, op. 

Laat zien dat de cycloïde door de 2 punten en die 

zijn top heeft in het bovenste punt hieraan voldoet. 

Schets de brachistochroon. 

Wat zou de brachistochroon zijn voor het analoge pro- 

bleem, waarbij echter de puntmassa een beginsnelheid 


Dd 
en krijgt ? 


Opmerking |. 
Een cycloïde is de baan die een punt op de omtrek van 














Instituut voor Theoretische Fysica - kh 
der Rijksuniversiteit Utrecht 


197k-1915 (vervolg 7) 


een rollend wiel beschrijft, zo'n curve heeft dus als 


Werkkollege Statistische Mechanica 





parametervoorstelling : 
y A(1 -— cos t) „5 
x= Alt - sin t) 


Welke parameters kunnen aangepast worden zo dat deze 


Li] 


cycloïde naast het punt (0,0) nog een ander gegeven 


punt (x, k En bevat ? 


Opmerking 2. 


De differentiaalvergelijking waaraan de brachistochroon 
voldoet kunt u in cartesische coördinaten terugbrengen 


tot een oplosbare eerste orde differentiaalvergelijking 
dy 

dx ä 

Ook kunt u de hele vergelijking vermenigvuldigen met Se e 


door te scheiden naarde variabelen y en p= 5 
waardoor er een direct integreerbare differentiaalverge- 


lijking ontstaat. 


Een elektrische puntlading met massa m bevindt zich in 
een homogeen élektrisch veld Ë en magnetisch veld Ë is 
beide niet van de tijd afhangend, en is bovendien door 
een ideale veer (veerconstente a) verbonden met de oor- 
sprong. 
Schrijf op de Lagrangevergelijking van de 2e soort en 
de Hamiltonvergelijking en los deze op. Beschouw ook 


in het bijzonder het geval a = 0. 


Een staafje ter lengte l en massa m, met verwaarloosbare 
doorsnede, is vrij draaibaar opgehangen in zijn zwaarte- 
punt, en bevindt zich in een tijdonafhankelijk magnetisch 
veld. De magnetische eigenschappen van het staafje worden 
voldoende beschreven door zijn dipoolmoment dat evenwij- 
dig staat aan de as van het staafje. (Alle hogere multi- 8 
poolmomenten zijn nul). 

Beschrijf de beweging van de 2e soort en de Hamilton- 

vergelijking. 

Herleid de vergelijkingen tof die van een 1-dimensio- 


nale oscillator. 


C En 


zl 
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MG (Viriaaltheorema ideaal gas). 


Beschouw een afgesloten vat{l met volume V, waarin N 
puntvormig gedachte molecuien met massa m,zonder wissel- 


wenking. Uitwendige krachten behalve die van de wand, 
crate miner tekracht) worden verwaarloosd. | 

Leid voor dit systeem, het zgn. wideale gas", uit het 
viriaaltheoremá, de wet van Boyle-Gay Lussac af. 
(Bedenk dat de gemiddelde kinetische energie per vrij- 


heidsgraad hier } k T is). 


Mi (ldeaal gas, Ensembletheorie, Kansverdelingen). 


Beschouw he tzelfde systeem als in opgave 10. 

Zij @ een deelvolume vanf) van willekeurige vorm, met 
volume v, =| oo | ‚ en zij n het aantal molekulen die zich 
op een gegeven moment in) bevinden. (n is een 'kans- 
grootheid" en wordt daarom onderstreept.) Wij kunnen 
niet n zelf berekenen, wel bijv. de kans dat n een be- 


paalde waarde n (n = 0,1 ‚2,...) heeft (Deze kans zullen 


wij Pn noemen) , mits wij een geschikt waarschijnlijkheids 


postulaat stellen, 
A) Kies zulk een geschikt postulaat 


B) Bereken Pp, en schets het verloop als funktie van n. 
(De "Binomiaalverdeling'). 


Bereken m = n ‚d.w.z. het gemiddelde aantal deel- 








tjes in , én de variantie van n , ot? ‚ uit 
or? =<(n - mi>, (TO = "standaarddeviatie"), 
(Aanw. : Bereken eerst de "genererende funktie'' 


oo 
2 Pa zh = p(z) en daaruit door differentiatie, m enor). 
nzo 


C) Beschouw de limiet Veco, Neo, maar Vo en ve y vast. 
Laat zien dat m , Gen P, elidig blijven, en dat in deze 
limiet P, nog steeds een (diskrete) kansverdeling 
voorstelt (d.w.z. dat nog steeds geldt dat 


Di B: * 1). Dit is de Poissonverdeling. Schets deze, 
n=0 


D) Welke limietovergang (verschillend van die onder C) be- 


schouwd) levert de Gaussverdeling op? 


E) Tracht voorbeelden te vinden van de limiet overgangen 
(of "benaderingen') die met C) en D} korresponderen, 
bijv. in het dagelijks leven, en/of in de ensembletheo- 


rie, 


AZ. Beschouw een Één-dimensionale harmonische oscillator. De 


kinetische en potentiële energie worden gegeven door 


| 2 
T-= Vel imuwexe. 


a) Bereken de met x gekonjungeerde impuls. 

b) Wat is de hamiltoniaan, Los de Hamiltonvergelijkingen 
op. | 

c) Teken de baan van het deeltje in de (twee-dimensionale) 
faseruiímte. 
Verifieer de stelling van Liouville. 

d) Geef de Maxwell-Boltzmann verdeling voor N deeltjes. 


Reken de toestandssom uit. 


e) Bereken de gemiddelde energie per deeltje. za 





an AS" indd an E75 ene „ed Lema ee ET 
es 5 een an eden 


hi tenet iernemorhede ichi 
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AS Beschouw een tweedimensionale harmonische oscillator met 


Omer ttondaan 
2 è 
H(p oP »x,y) = ae + mw (x + yy). 


Volgens het equipartitietheorema is de gemiddelde energie 
gelijk aan 2kT. Ga over op poolcoördinaten: 

Xx == r Cos & 

y= r sin 6 « 

Bereken de net r en ’ gekonjungeerde impulsen Pp, en P, en 
laat zien dat de Hamiltoniaan gegeven wordt Rade. 


2 
H(p P,ereó) = ES Ë hi ze | zino 
” 
Er zijn nu twee kinetische-energie-termen en één kwadra- 
tische potentiëfe-energieterm. 
Verifieer door direkte berekening in pooifkoördinaten, dat 
de gemiddelde energie gelijk is aan 2kT. 
Wat betekent dit resultaat in verband met het equipar- 


titietheorema ? 


WS ve beschouwen een ideaal gas in een zwaartekrachtveld 
V(z) = mgz. | 
m is de massa van de deeltjes. We bekijken een zuil met 
Os x <«L 
Os ysL 5 
0 << zg e 
a) Leid de barometrische hoogteverdeling af door de 
krachten op een dun laagje van het gas te berekenen. 
b) Geef de Maxwell-Boltzüann-verdeling voor N deeltjes. 
Bereken de snelheidsverdeling op de hoogte z, en de 88e Zee 
middelde snelheid. ek 
Wat is de gemiddelde energie per deeltje op de hoogte 
Z, en wat is de gemiddelde energie per deeltje? 
Wat is de gemiddelde hoogte van de deelt jes? 








c) Ons interesseert alleen de beweging in de z-richting. 
Voor het gemak bekijken we nu het één-dimensionale 
ideale gas, waarbij de plaats van een deeltje gegeven 
is door de koördinaat 3 de impuls door p. De poten- 
tiële energie is dezelfde als boven. Los de bevegings- 
vergelijkingen op bij de beginvoorwaarde (Poszo): 

Geef de vergelijking voor de baan van een deeltje in de 
(tweedimensionale) faseruimte. 
Teken een paar mogelijke banen in de faseruimte. 


15. a) Wat is (volgens Maxwell) de kans dat een deeltje een 
snelheid heeft die ligt tussen Ven Vv + dv (hiermee 
wordt bedoeld dat de x-component ligt tussen Fe en v‚ + 
dv, » et C Ns 


b) Wat is de kans dat een deeltje een snelheid heeft 


waarvan-de-tengte Ligt-tussen-v + dv. » as 


c) Wat is de kans dat de energie van een deeltje ligt tus- 


sen E en E + d E , (als de energie alleen uit 


kinetische energie bestat)? z * 


d) De inhoud van een bol in een n dimensionale ruimte is 
ze 
FUTastT RP, waarin R de straal van de bol is. 

Bereken nu de kans dat de totale energie van v onaf 
hankelijke deeltjes ligt tussen E en E + d E , waarbij 
veer de potentiële energie verwaarloosd. mag worden. 
Bereken de gemiddelde energie, en de spreiding in de 
energie voor v onafhankelijke vrije deeltjes. 


e) Gegeven 1 mm lucht bij kamertemperatuur, wat is dan 
volgens deze berekening de relatieve spreiding in de 
energie van dat systeem? (De temperatuur wordt exact 
bekend verondersteld.) 
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f) De kansen berekend onder c) en d) kunnen respektievelijk 
geschreven worden als p (E)4E en p,(E)dE, 
1 
Bereken voor algemene waarde van v(v=1,2,3,...): 
1) De energie E ax waar voor p,(£) een maximum heeft. 
2) De gemiddelde energie E . 
gem 
Laat zien dat voor grote waarden van v het verschil 


… E | 
gem max 


hoogstens van dezelfde orde van grootte is nls de stan- 

daard-deviatie oc van de verdeling (Co, (H)ar) van de energie 

E. (Dwz. bewijs dat |E „- E begrensd is als v>e), 
0 gem max 


g) Bewijs dat 


+ oo _ 2 N 
fe ax ax e= Á (als a>Q) | (1) 


me OO 
Leid hieruit af dat RE 
r(3) = Am 


Verhef nu beide leden van (1) tot de ate macht. (Schrijf 
het linkerlid als een n-voudige integraal). Herschrijf nu 
het linkerlid in termen van poolkoördinaten, en bewijs 
uit de zo gevonden relatie de onder d) vermelde formule 


voor de inhoud van een n-dimensionale bol. 


BAE nog eens door, de gronábeginselen van de Thermodynamika, zoals 
bijvoorbeeld beschreven in : Becker, "Theorie der Wärme "‚JÔ 1-12, 


N 13 a), 619 en Ó20, of in Fermi, "Thermodynamics", 9) 1-14 en9Q16-18. 


17. Dipool-dipool wisselwerking tussen atomen. 
We beschouwen twee atomen (ongeladen), ieder bestaande 
uit een zware kern en één elektron. De elektronen zijn 
harmonisch gebonden aan hun kern (frequentie w)- De af — 
standsvektoren van de elektronen tot hun kern zijn resp. 
P, en En Beide kernen zijn gefixeerd op een onderlinge 
afstand R (langs de z-as bijvoorbeeld), en vormen een 


twee-atomig molekuul (b.v. H‚). 





ahd 
ahd 


5 ' aA 





Als we de atomen opvatten als dipolen met dipoolmomenten 
u, had er, en uz = er, ‚ dan kunnen we de wisselwerking 
tussen de atomen beschrijven als die tussen twee dipolen. 
De potentiële energie van een dipool met dipoolmoment u, 
in het veld Ë, van een andere dipool is gelijk aan -u,Ën- 
Het veld van een dipool met dipoolmoment A wordt gegeven 


door 





u Rv 
É,  … Â + 3 5 R 5 
R R 


waarbij R de afstand tot die dipool is. 


a) Schrijf de volledige Lagrangiaan op voor dit systeem. 


b) Voer de volgende transformatie uit: 








/e 
en bereken de met q, en do toegevoegde impulsen B, en 
=>} 
Bs « 

c) Geef de Hamiltoniaan H(P, Po sd, vdo) en de Maxwell- | 
Boltzmann verdeling voor N van deze twvee-atomige mole- 


kulen. 
d) Reken de toestandssom en de vrije energie uit. Bepaal 
vervolgens: de gemiddelde energie per molekuul en de pe 


soortelijke warmte Cc, per molekuul. 

Reken de gemiddelde kracht uit tussen twee atomen van 
een molekuul. Met welke macht van R gaat deze kracht? 
Is zij aantrekkend of afstotend? Deze kracht is de 
Van der Waals kracht. 
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18% Een grammolecuul van een twee-atomig ideaal gas (tempera- 
tuur 18°C) zet adiabetisch uit tot 1,35 maal zijn oor= 
spronkelijke volume. Wat is de eindtemperatuur en de door 


het gas verrichte arbeid? {R = 2,0 cal/graad, 1 cal = 
b,2,107 erg). 


197 Een grammolecuul van een ideaal een-atomig gas verricht 
een Carnot-cyklus tussen de temperaturen ho0®K en 300°K, 
Bij de isothermische toestande verandering bij de hoogste 
temperatuur is het beginvolume 1 liter en het eindvolume 
5 liter. Bereken de hoeveelheid arbeid verricht gedurende 
een cyklus, en de hoeveelheden warmte die door de twee 
warmtereservoirs gedurende een cyklus worden opgenomen 


resp. afgestaan. 


Be. Bereken het minimale vermogen dat een koelkast verbruikt, 
indien de inwendige temperatuur -— 3% is, die van de om- 
geving art, en indien er een hoeveelheid warmte van 5 


kcal per minuut naar binnen lekt. (1 cal = hb Joule). 


AT Ben grammolecuul van een twee-atomig ideaal gas ondergaat 
een toestandsverandering, welke in het (p‚V)-diagram door 
een recht lijnstuk wordt beschreven. In de begin- en 
eindtoestand zijn de temperaturen resp. in en Te en de 
volumina resp. V_ en V . Bereken de door dit systeem ver- 

‚richte arbeid en geabsorbeerde warmte, uitgedrukt in 


Ts V‚ (id = 1,2) en R. | 


ee 


22, Beschouw het systeem van opgave 6, waarbij echter de ke- 
teh van massa's en veren wordt uitgebreid tot n massa's 4 
en n + 1 veren. Bepaal weer voor kleine trillingen de 
eigenfrequenties van dit systeem, en geef de verdeling 
van die eigenfrequenties in de limiet n*e , Doe dit op ‘ 


twee manieren: 











(vervolg opg. 22) 


I. 


II. 


„de (gegeneraliseerde) orthogOnaliteitsrelaties. 


(Via constructie van de eigenvectoren). 
Probeer als uitwijking van het j-de deeltje: 
u(t) = a, eiet , 
d J 
Vind m.b.v. de bewegingsvergelijking een verband tussen 


a etc. 


dT gert Mien, 
Substitueer a, = AÌ R 

en bepaal hieruit en uit de randvoorwaarden dat de muren 
vast zijn, w en À 

Controleer of het aantal vrijheidsgraden gelijk is aan 
het aantal zo gevonden eigentrillingswijzen. Is het ge- 
vonden stelsel lineair onafhankelijk? Volledig? Verifieer 


(Via de eigenwaardevergelijking det (V-u°T) = 0). 
Noem deze determinant fa ‚, Vind een relatie tussen Ld n 
bn en B IE Herleid deze, door geschikte keuze van para= 


meters y en x,‚ tot de volgende recursierelatie voor de 


grootheden 8, =Y A 


-2xe, + En =O: (1) 


Laat zien dat hieraan voldoen de Tchebitchef polynomen 


En+1 


TG), maar ook de verwante polynomen Ss), gedefinieerd 


„door 


Tr) = cos(n arc cos x) 


Ax sx) = sin (n are cos x), 


en dat de meest algemene oplossing van (1) geschreven kan 


worden als 
Ex) = a(x) TG) + b(x) Ss (x) 


voor zekere a(x) en b(x). Bepaal a(x) en b(x) uit de Ì 


beginvoorwaarden 8, (x) E «. en &,(x) u ee 


hl 
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23. 


25, 


(Vervolg). Ter vermijding van misverstanden in verband met 

vraagstuk 17, i.h.b. gedeelte dà), nog de volgende opmer- 

kingen: 

le) Bij molekuulvorming (bijv. EH) spelen Van der Waals- 
krachten een te verwaarlozen rol. (Welke zijn wel be- 
langrijk?) 

2e) Het aangegeven model is wèl goed bruikbaar om de Van 
der Wsalsikentht uit te rekenen tussen twee deeltjes 
die (in chemische zin) niet tot één molekuul behoren 
(Bijv. tussen twee A atomen, of twee N -molekulen), 

3e) In dit vraagstuk worden, uitsluitend met het oog op de 
berekening, de twee wisselwerkende atomen even als één 


“molekuul" opgevat. 


Beschouw, van het systeem van vraagstuk 22 een groot aantal, 
onderling niet wisselwerkende replica's die in evenwicht. 
zijn met een warmtebad bij temperatuur T, 
a) Bereken de soortelijke warmte per systeem. 
b) Bereken hieruit, in de limiet n+e de soortelijke warmte 
per deeltje. 
ce) Vergelijk de zo gevonden uitkomsten (wat betreft 
frequentiespektrum en soortelijke warmte per deeltje) 
met de aannamen hierover resp. de resultaten hiervoor, 
in de Debije- en in de Einsteintheorie van de vaste stof 


in 1 dimensie, 


Bereken de soortelijke warmte van het tweedimensionale 


Debijemodel. 


Veronderstel dat, in plaats van een dispersierelatie van de 
vorm w = slal (waarin w de frequentie, a de golfvektor en s 
een konstante is), welke geldig is voor een vaste stof vol- 
gens het model van Debije, er een dispersierelatie geldt 


van de vorm: 


w=clal® 


waarin c en n konstanten zijn. 








(vervolg opg. 25) . ie me Th 
‘ Bereken het gedrag, bij lage temperatuur, van de soorte- 
lijke warmte van een vaste stof met zulk een dispersie- 


relatie, 


‚ De energie van éen tweedimensionale harmonische oscillator, 
als funktie van twee kwantumgetallen Et, en Er wordt ge- 


geven door: 
€ __= hw (n + 3) + hw (n + }) 
1 1 2 2 


n ‚n = 0,1,2,.s« 
2 


a) Bereken de kwantummechanische toestandssom. 

b) Bereken, voor een groot aantal identieke exemplaren van 
deze oscillator, zonder onderlinge wisselwerking maar in 
evenwicht met een warmtebad bij temperatuur T, de vol- 
gende grootheder: 

1) De vrije en de gemiddelde energie per oscillator. 
2) De soortelijke warmte per oscillator. 

c) Veronderstel nu dat w <<w , en teken c als funktie van T, 
Vergelijk dit met het Eimateke resultaat. Wat betekent Ì 
dit voor het equipartitietheorema? 

d) Bereken de toestandssom voor de isotrópe oscillator 


(wo == w ):- 
1 2 


le) door in het resultaat van a) u, = w te nemen; 
2 
2e) uit het energiespektrum, door expliciet de ontaar- 
dingsgraad van elk energieniveau in rekening te 


brengen. 
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/ 
27. Beschouw een atoom in de grondtoestand, met spin 3 . 
PA 


id Laat hierop een magnetisch veld B langs de z-as werken, 
(Zwak genoeg, zodat hoger aangeslagen toestanden van het 
atoom buiten beschouwing gelaten kunnen worden). 

a) Bereken de kwantummechanische toestandssom. 

b) Bereken, voor een groot aantal identieke exemplaren van 
dit systeem, zonder onderlinge wisselwerking maar in 
evenwicht met een warmtebad met temperatuur T (alles 
per atoom): 
le) Het gemiddelde magnetische moment langs de z-as, me 
2e) De inwendige energie u. 
3e) De soortelijke warmte c. 
he) De entropie s. 

c) Schets het verloop van deze vier grootheden als funktie 
van T, Leid asymptotische uitdrukkingen af‚ geldig in de 
limietgevallen T+o en T»e, en ga na of deze resultaten 


ook op andere gronden te verwachten zouden zijn. 


„28. (Thermo) 


a) Bewijs dat de chemische potentiaal u = 5 gegeven wordt 
door % 
ADD) 
T‚v 


b) Bereken voor een ideaal gas u als funktie van p en T. 
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J De ‘toestandvergelijking van Van der Waals, welke bij benadering geldt 


voor een klassiek niet-ideaal gas wordt meestal als volgt geschreven: 
(p + 27 (V-b) = RT 

(Deze formule geldt voor | grammolekuul van het gas). 

Hierin zijn p de druk, V het volume, R de gaskonstante (= N‚k „ke 

konstante van Boltzmann, N, = getal van Avogadro), T de temperatuur, 
“en a en b konstanten die karakteristiek zijn voor het beschouwde soort 

gasmolekulen. 

Schrijf deze formule in de vorm van de zgn. 'viriaalreeks' (dwz. reeks- 


ontwikkeling van p in machten van p ‚ het aantal deeltjes per volume- 


eenheid in het systeem) en bereken de eerste drie viriaalkoefficienten. 


30. We bekijken een (klassiek) ideaal gas (N molekulen) in een zwaarte” 
krachtveld gegeven door V(x,y,z) = m g z. Het gas bevindt zich in een 


oneindig.hoge ruif, gegeven door 


0 <x<ztg a 
0 < y<b 


0 < ez 





a) Geef de Maxwell-Boltzmann verdeling. 

b) Bereken hiermee de (plaatsafhankelijke) deuk p(z). 

c) Bereken het (moment vanhet) koppel dat werkt op het gearceerde vlak. 

d) Bereken de toestandsintegraal.  … it \ 

e) Bepaal de vrije en gemiddelde energie per deeltje,,en de soortelijke | 
warmte bij konstant volume, cy: 

£) Bereken wederom het (moment van het) koppel dat werkt op het gearceerde 


vlak, maar nu uit het resultaat van d). 


€ % í at NN 








EE 


31. Bewijs, voor een kanhoniek ensemble, van een Byuteem van 
N deeltjes met linten Helle, b), dat 


<A, > = _ dr 


da, 
28F 
z t ede 
U z= <H> Tan 
<a)? E uk. 
kr 


da; Es 


a 1 F 
elan han te” ek” va: ze | 


(Hierin is A, = === , AH « H … eH>, AA, = A, - <A, > . 


<AH AA, > = 


Dj 


De grootheid A, is de tgn. “inwendige parameter” behorende 
bij de waitwendige parameter” a, 
A, en a, vormen samen “een parameterkoppel"). 


„à eN 
Merk op dat Ta; <> id Ta, He : 


Bereken alle bovenstaande grootheden voor het systeem be- 
schouwd in vraagstuk 30, Beschouw hierbij als uitwendige 
parafeters : a , Db , & 

Wat stellen de bijbehorende inwendige parameters voor fysische 
grootheden voor? 


32 , In een keten van N molekulen kan elk molekuul zich op 
twee manieren oriënteren: Als het in de toestand a is, 
heeft het een lengte a en een energie B 8 als het in de 
toestand B is, heeft het een lengte b en een energie E 
Op de uiteinden van de keten werkt een spânning S. De 
lengte L van de keten is gelijk aan de som van de lengtes 
van de samenstellende molekulen. Leid een relatie af tus- 
sen S en de gemiddelde waarde van L bij een temperatuur T, 


B 





2} 
Bereken de relatieve fluctuatie van Li sale 
in de limiet Nee , 
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35. Beschouw een systeem met drie energieniveaux 0 , e en 2e waarvan het 
middelste tweevoudig oataard is en de andere twee niet ontaard. 


Bereken de soortelijke warmte van dit systeem. 


Á. Klassificeer de volgende deeltjes (fermionen, bosonen). 


a-deeltjes; 3e 8 bne s H 53 H, : e“ (positon) ; 
f 


‚+, Jest 
Li =ion ; Li =ion . 


35, Gegeven: 3 onafhankelijke deeltjes in 3 {l-deeltjes-toestanden. 
Verwaarloos de energieverschillen tussen die toestanden. 
Bereken de kans, in de drie statistieken, op de volgende mogelijke 
gebeurtenissen. | 
a) De drie deeltjes zitten alle in dezelfde toestand; 
b) Twee zitten in eenzelfde toestand, één in een andere; 
c) Elke toestand is bezet. 
Welke statistiek zoudt u het gezelligst noemen, en welke het onge- 
zelligst? 
„3 Bereken hoe het aantal bosonen in de prondtoestand, van een 
ideaal Bose-Einsteingas met een temperatuur beneden de kon- 
p densatietemperatuur, van de temveratuur afhangt. 
‚37, Bewijs dat een twee- en een eendimensionaal gas van niet- 


wisselwerkende Bose-deelt jes geen Bose-Einsteinkondensatie 


vertoont. 


Bereken de struktuurfunktie u(E) = f6(H, - E)dr van een 
klassiek ideaal gas van N deeltjes (zonder inwendige 
vrijheidsgraden) in een volume V. 

(Bereken eerst de kanonieke partitiefunktie-voor dit 
systeem, en maak gebruik van het algemene verband tussen 
de struktuurfunktie en de kanonieke partitiefunktie), 


Pe 


Pe Beschouw een systeem van. N onafhankelijke quantummecha- 
nische rotators met. Boltzmannstatistiek, Het traagheids- 
moment van elk deeltje is JT. 





2 
Ee 5 eel) ‚ ontaarding g(l) = 21+1 . 
è 
Ö _ A 
Voer in 8 = BIKk * 


kies j 





ho. 


De toestandssom is gelijk aan 


- 19 - 


Sereken de toestandssom voor \ -2 
a) T ea e f 5 ó hit 
an WE 
b) T >> 6 NE 
kt 7 ee 3 
Bereken in beide gevallen de soortelijke warmte C(T). ok 


Maak opgave 39/nogmaals maar nu voor zowel Fermi- als 
Bose-statistiek, en met N=2. (Gebruik niet het grand 
ensemble). 

Bereken ook de groot-canonieke toestandssommen, van beide 
statistieken en bepnal ook daaruit weer het vorige resul- 


taat voor N=2, 


Twee volumina I en II bevatten samen N deeltjes en zijn 
in thermisch kontakt met een warmtebad van temperatuur 
T. Er bestaat een of ander mechanisme waardoor de twee 
volumina deeltjes kunnen uitvisselen:;N = N, + Ni: = 
constant. 

Wat is de voorwaarde voor de meest waarschijnlijke ver- 
deling? 

Wat is de waarschijnlijkheid dat N, deeltjes in I en Ni: 
deeltjes in II zitten, voor het geval van een ideaal gas? 


De Ising keten. 

We beschouwen een keten van N vastgeprikte spins. Iedere 
spin, gekarakteriseerd door de variabele 8, » kan de waar- 
den s, =+ 1 (spin up) sn 8, 1 (spin down) aannemen. 
(i = 1 „2,ese,N). De énergie van een configuratie (s, } 
wordt gegeven door 


KT iks J sisser 


del De 


= L exp {B J d 8, 5; K 5 B 


waarbij Kerohlieùta wordt over alle configuraties ts, } « 
Bereken de toestandssom, de gemiddelde energie en 


de soortelijke warmte. 
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h3, Beschouw twee volumina v, en Vo ‚ die verbonden zijn door een buis. 
Het geheel bevat N onafhankelijke deeltjes. Tussen v, en V‚ heerst een 
potentiaalverschil zodanig dat de potentiële energie van een deeltje 
in V‚ gelijk is aan nul, en de potentiële energie van een deeltje in 


v, gelijk is aan W. Het geheel bevindt zich in evenwicht. 


a. Bereken de toestandssom, de vrije energie en de chemische potentiaal 
van dit systeem. 

b. Bereken de druk in v, en de druk in Vo: In hoeverre geldt de wet van 
Boyle-Gay-Lussac? 

c. Bereken de gemiddelde aantallen deeltjes N, en N, in resp. v, en Vo. 

d. We sluiten de verbinding tussen v, en Vo: Bereken de chemische po- 


tentialen hj en Ho in v, en Vo en laat zien dat Ho + W= Ho 


HC ve hebben N deeltjes verdeeld over M roosterpunten (M>N). Ieder deeltje 
moet zich in een roosterpunt bevinden, en ieder roosterpunt kan maxi= 
maal één deeltje bevatten. De toestandssom voor dit systeem is 

| ZM) = Kew 
waarbij W het aantal mogelijke manieren is waarop N deeltjes M rooster- 
punten kunnen bezetten. (Merk op, dat de waarschijnlijkheid voor iedere 
bezetting even groot genomen is). 
a. Bereken W. 
b. Bereken de vrije energie F „(Ms B). 
ce. Neem N en M zeer groot hat. Bereken nu Ky door te eisen dat de 
vrije energie een extensieve (= additieve) grootheid is. 
d, Bepaal nu de chemische potentiaal u. 
e. We identificeren M met het volume. Bereken de druk. 
f‚ Bereken nu de Gibbs-potentiaal G(N,p,T) en verifieer G= u N. 


g. Bereken de entropie en verifieer S = k log W. 





L5. 


L6, 


_ 21 -— 


Op het kollege is de grand-potentiaal berekend voor het quantummecha- 
nische systeem van niet wisselwerkende deeltjes (zonder spin) met 


Maxwell-Boltzmann statistiek: 


RE =Be, 
ala,8,v) =e Je 
Vv 


… . Le KJ 
waarbij gesommeerd wordt over alle energiewaarden van één deeltje; Ee, 


Wij willen deze uitdrukking nu schrijven 


ala,B,v) = fe FD(vae 
le} 


welke uitdrukking moet gelden in de limiet dat het volume V naar on- 
eindig gaat. 

Bereken de toestandsdichtheid D(e). Wat is D(e) voor deeltjes met 
spin s ? 

Leid ook de overeenkomstige uitdrukkingen af voor q(a,8,V) in het ge- 


val van Bose-en van Fermistatistiek. 


Beschouw N deeltjes. Ieder deeltje kan Òf harmonisch gebonden zijn (met 


„frequentie w en een extra bindingsenergie =X) of zich vrij bewegen in 


een volume V, 
We hebben dus een systeem van een (veranderlijk) aantal driedimensio- 


«nale harmonische oscillatoren met frequentie w en een veranderlijk aan- 


tal vrije deeltjes. De som van het aantal vrije deeltjes en het aantal 
oscillatoren is gelijk aan N. Het geheel is in evenwicht. 


a. Bereken de kanonieke toestandssom Z(B8,V). Welke evenredigheidsfak- 


tor moet u toevoegen teneinde de bride energie een extensieve 
grootheid te maken. 

b. Bereken de groot-kanonieke toestandssom en de grand-potentiaal. 

c. We beschouwen dit systeem als een model voor een damp en een vaste 
stof in evenwicht. De damp wordt beschreven als een ideaal gas, de 
vaste stof wordt beschreven als een aantal identieke drie-dimensio- 
nale harmonische oscillatoren zonder wisselwerking. 

(Einstein model). Bereken het gemiddeld aantal deeltjes in de damp 


<N> en in de vaste stof <N> ‚ bij gegeven temperatuur en chemische 


D 
potentiaal. 
d. Schets <N> als funktie van 8 = EE: Wat is het gedrag voor 


/<N> 
T>0 ? Ny 


Wat is het gedrag voor Te ? Kunt u een verklaring geven voor dit 


merkwaardige resultaat? 
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(vervolg opgave h6) 
e‚ Bereken de druk voor een vast totaal aantal deeltjes <N> + <N> 


als funktie van T. 


h7. Beschouw het systeem van opgave h6 en bereken hiervoor de quantummecha- 
nische grand-potentiaal. Neem de deeltjes onderscheidbaar. De bindings- 
energie van de deeltjes in de vaste stof is groter dan de nulpunts- 
energie (per deeltje): 


X > ERS w 


Los weer de onderdelen c., d. en e. op. 


48, Beschouw een systeem van niet-wisselwerkende deeltjes waarin het maxi- 
maal toegestane aantal deeltjes per energie-niveau gelijk is aan d. De 
deeltjes zijn ononderscheidbaar. 

Bereken de grand-potentiaal voor dit systeem en de gemiddelde bezetting 
van de niveaus. Kontroleer de uitkomsten voor d = 1 en d= @, 
Opm, Deze statistiek die niet in de natuur voorkomt (afgezien van d = 1 


en d = ©) wordt wel "intermediaire statistiek genoemd. 


19, De definitie van de g-funktie van Riemann is: 


oo 


cls) = } L (Re s>1). 
1=1 1 
We geven: 
5 e | | 
c(5) = 1,3k2 | 
2e | 
E(5) = 2,612 . 


8,(a) = 


Bewijs met het bovenstaande de eerste drie termen uit de volgende asymp- É 
totische ontwikkeling: | 


85 jole) = 2,363 (-a) 2/2 + 1,342 + 2,612,a 
_= 0,730.a°+ (a°) EE: 











“0e 
50. Beschouw een ideaal Bose-gas. Bewijs dat in het overgangspunt Te de 


afgeleide van bs naar T een sprong maakt die wordt gegeven door 


Lm 5) … (B Ze) … 3266 
ark NT, ar Nk MT, Te 


51. Voor een kristal bestaande uit N atomen, elk met z elektronen, kan men 


definiëren: 
kr = de straal van de Fermibol (in de theorie van de vrije elek- 
tronen) 
a = de maximale waarde van het golfgetal van het trillingsspektrum 
(in de theorie van Debije voor de kristaltrillingen. Longitu- 
dinale en transversale voortplantingsmethoden gelijk genomen). 


Vind het verband tussen k 


pen ps uitgedrukt in z. 


52, (Het kwadratische Stark effekt) 

Beschouw een halter (zie. opgave 1) bestaande uit twee gelijke massapun- 
ten m, op onderlinge afstand a, met ladingen te en -e , in een elektrisch 
veld F, in de z-richting. 
(Traagheidsmoment I= } ma° ‚ elektrisch dipoolmoment p = ae). De halter 
is vrij draaibaar maar vastgeprikt in zijn zwaartepunt. 
Bij een quantummechanische behandeling van dit systeem vindt men, in 2e 
orde storingsrekening, de volgende uitdrukking voor de gestoorde ener- 
gieën van dit systeem: 

Bin = zien + 3 bn r‚° + or,» 
waarin de quantumgetallen j en m beide geheel zijn, en voldoen aan 

j>o,-jSmSj. | 
De coefficiënten Ber kunnen ook berekend worden, en zij blijken te vol- 
doen aan 

Je Den =oals j #0 


terwijl boo gegeven is door 


Beschouw nu een systeem van N van zulke halters, zonder onderlinge wis= 
selwerking, met M.B. statistiek. Bereken in dezelfde benadering (aus 

tot in 2e orde in F‚) de gemiddelde waarde <P> van de z-component van 
het totale elektrische dipoolmoment Ë ‚ en hieruit de elektrische sus- 


ceptibiliteit Bor 
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(vervolg opgave 52) 


"Di 


Geef ook alle andere componenten van de tensor X 
Schets het verloop van Ke 0) als funktie van T, Beschouw in het bij- 
zonder de limieten Tso en ‘le en geef in deze Limieten benaderings- 


formules voor Xez 


Bereken de soortelijke warmte in 2e orde in F,. 


53. 


Beschouw het volgende model voor een systeem van polymeren. Een poly- 
meer ter grootte n bestaat uit n atomen A met massa ue Een polymeer 
bestaande uit n atomen heeft een bindingsenergie Xn De polymeren be- 
wegen verder als vrije deeltjes met massa m, Sn} door het volume V. 
Als er N, atomen A zijn, bewijs dan dat het aantal mogelijkheden om 


a) a 
deze in N, polymeren ter lengte 1, N, polymeren ter lengte 2,..., 
N, polymeren ter lengte n‚.…. te verdelen gelijk is aan 
ee De bs 
Ï Nn Î a ie 
Et 
n=1 nej 
Merk op, dat 
co 
J n.N =N, 
nst Ans An ha ER 
b) Bewijs dat de, kanonieke toestandssom gelijk is aan 
as if der 3/2 N 
: L 3 Ë k 
Ket re Er, 
A nei “n° ‘n!h 


waarbij de som gaat over alle verdelingen {N_} met 


A 
ami n A 


Cc) Bereken de groot-kanonieke toestandssom en leid hieruit af de druk 


hd 


en het gemiddelde aantal atomen A bij een chemische pven tian Hye 
d) Wat is het gemiddelde aantal polymeren <N > met grootte n? Bereken 
ook de gemiddelde fluctuatië in dit aantal.” 
e) Wat is het gemiddelde totale aäntal polymeren 


<i> 
n= 
f) Bewijs, dat En N. Had u dit ook verwacht? 
g) Wat is de gemiddelde grootte van de polymeren? 


Merk op dat dit model ook gebruikt kan worden als ruwe benadering voor 
een damp. De polymeren ter grootte n corresponderen dan met een druppel 
bestaande uit n atomen. : 
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Vrije elektronen bezitten een permanent magnetisch moment gelijk aan 
u. e Ha . Kd 


waarin ug (= 1 Bohrmagneton) en 10720 erg/gauss, Br soin van het 
elektron. In een uitwendig magnetisch veld worden de magnetische momenten 
gericht; de bijdrage van de epin van de vrije elektronen geeft het Puaulí- 
spin paramagnetisme. En 
Laat B een homogeen uitwendig magneetveld in de z-richting zijn. De spintoe- 
stand van het elektron kan beschreven worden door de z-komponent S,. Het mag- 
neetveld geeft een energieverschuiving 


nev -RiE- agee Sujke 


korresponderende met de apintoestanden s, =t} . De energteniveaux van een 
vrij elektron met spin zijn 

22 

era = om “vat 
Hierbij is de wisselwerking tussen het magnetisch veld en de baanbeweging ver- 
waarloosd. (Deze wisselwerking geeft aanleiding tot het Landau-diamagnetisme). 
a. Bewijs dat de magnetische susceptibiliteit van N vrije elektronen in een 
volume V voor zwakke aanbad gegeven wordt door 


x= 2u5 [ de D' (e) CEP 
ééndeeltjes … 
D(e) ís de Vniveaudichtheid van een gas van vrije fermionen zonder spin. 

b. Geef een benaderingsformule van X geldig voor hoge temperaturen. 

c.‚ Bereken X voor T =O, en vind bovendien de eerste korrektieterm. (Maak daar- 
toe gebruik van de formule van Sommerfeld. Zie bijv. ter Haar, p. 152 of 
Becker, 2° druk, p. 171) 

Geef numerteke schattingen voor de fermi-energie van 

a. elektronen in een metaal 

b. nucieonen ín een: zware kern 

Ce SHe-atomen in vloeibaar. Sje 

(volume per atoom: h6,2 AS). 

Bereken de overeenkomstige grootheid voor bre. 

Breng deze in verband met fysische grootheden. 


Pl 


„5 


Bewijs dat een ideaal Fermigas voldoet aan de betrekking 
è 


waarin U de: nine energie van het gas is. Hoe groot is U bij het absolute 
nulpunt? 
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57. (Viriaaltheorema niet=ideaal gas). Hetzelfde systeem als 
in opgave “10. 5 ‚ behalve dat nu de (puntvormig of bol- 
vormig gedachte) deeltjes wèl krachten op elkaar uit 

oefenen. Potentiële energie van de wisselwerking 


Ee « + 
ï ellr, - rl) é 
tei sjen 
Hierin zijn P, - Fn de coordinaten van de middelpunten 
van de gasmoleculen, en é(r) de "paar-potentiaal", is 


KJ } 
een gegeven funktie van r. 


gf 
Voor de meeste gassen ziet 


é(r) er ongeveer zo uit: 





<afstotiag a entrakt ende det 
A) Splits de totale viriaal V in twee gedeelten V = jn + 
Vv; : Vo te&.v. de krachten v/d wand en Vv; tege ve ie 
onderlinge wisselwerking. 
Voer in de funkt ie s(r, ‚ F‚) ("paardistributiefunktie") 
als volgt: -. - 
Gegeven twee volume-elementjes är, en ar, 
ter plaatse van É, resp. E. ‚ dan is de 
kans dat in beide volume-elementjes zich 





een (middelpunt van een) gasmolecule bevindt evenre- 


dig aan odr, en aan odf, (os Fn). ge kde 
faktor is e(r, ni Fo ). Dus: kas = 0 telE, - Far, ar. 
Tracht ER te maken dat g de diene ele 
schappen heeft: . 

- 1) Als r. en F. ver van de wand, dan e(r, mn Fo) > 
2) Als |F, - False , dan slr, - Folder . 
3) Als é(r) = 0 dan a(r, é Ps ) = 1 (voor alle r, s F,). 

« 


kb) Als p+0 , den a(r) + e” el (Boltzsmannfaktor voor 
2 deeltjes) 





LAN 





/ 
58e 


= 27 — 


B) Druk de gemiddelde inwendige viriaar v, uit in de 
funkties élr) en glr), en leid uit het rbrlmnteneerdds 
af: 

2 Den 
en ep Ber Íre'(r)gtrdar (1) 
c) Bewijs hi ier vhn 


Er =o+ 2 jat y(r) (esr « el) (ze 


vaarin 


yr) = e en ar) 





(Gas van harde bollen). SO 
Systeem hetzelfde als in opk. 57 behalve dat wij voor 


è(r) de keuze maken: 
Lw) als O0 < r € ó 
lr) =pn (r)e egen 
ne als r > og 


waarin oe EN: ve. harde bol. 


dervank in gedachte de bollen door massapunten, elk om- 
geven door een "werkingssfeer" van straal Je 


Geen enkel massapunt mag zich bevinden binnen ue wer. 


kingssfeer van een ander massapunt. 


Noem P, de (gemiddelde) druk op de werkingsifeer’ van een 


willekeurig gekozen deeltje, uitgeoefend door de andere 


massapunten. 


A) Leid af, uit het viriaaltheorema, dat 


B - + by at 7 aks (1) 
k Ë Tof kT ĳ 
N 


(eigenvolume v.e. harde bol). 


B) Men kan aantonen dat voor harde bollen g(x.) een dis- 
continu verloop heeft in de Buurt van r = og , maar 


dat y(r) daar continu is: a 
1 





A 
“en 0740) 
pf ’ 
Dit aannemende, bewijs dan uit (2) van opg. 12 dat 


voor harde bollen: 


EE =p + hv oP 2 e(o+0) (2) 
waarin e(o+0) = Lim gote) ("Limiet van rechts"), 
e+ 


| 
“4 
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er ves (1) en (2) zien wij dat | | 

zr“ e &lo+0) OR 
Het rechterlid is de kans per volume-eenheid om een 
massapunt aan te treffen vlak buiten de werkingssfeer 
van een gegeven massapunt. 
Tracht de relatie (3) ook direkt af te leiden uit wat 
u weet over de Maxvell-Boltzmann snelheïdsverdeling 


van molekulen, 


D) Maak nu de veronderstelling 


P, … P 


en vind zo, uit (1), de benaderende uitdrukking voor de 
druk van een harde bollengass: 


Be 


T° lebvp à 


Ll 


59, Beschouw een systeem van harde bollen met diameter oc in 
à dimensies (d = 1, 2, 3). 
a) Bereken de 2e en 3e viriaalcoefficient (Lastig! ). 
b) Bereken de funktie g,‚(r) welke voorkomt in de dicht 
heidsontwikkeling van de paardistributiefunktie glr): 


„lr 
KT 


g(r) =e {1 + p e,(r) + p°e,(r) + eee} 


c) Bereken uit het resultaat van (b) op twee manieren de 
tweede viriaalcoëfficiënt . (Gebruik makend van het 
viriaaltheorema en van de uitdrukking voorde iso- 
thermische compressibiliteit),. 

à) Leid uit het resultaat van (a), via deze laatste zgn 
pressibiliteitsrelatie" een conditie af waaraan de funktie 


Ep 


Vcom= 


galr) moet voldoen. 








a 
/ 
| 








be 
66.* ï Een atoom met twee niveaux (encergieversch t 
ek e= Rw) bevindt zich in een ruinte met 
zwarte straling van tenperatuur T, De kans 
Ò per tijdseenheid op absorntie Zid wj‚g » 3 
en die op gestimuleerde plus spontane emisste A ‚ Neem aan dat geldt 
in „Be ‚a i ä 
Fo‚1 e ‚ ú Woo „ . ie 
A. Beschouw N, atomen waarvan N in het aangeslagen 
niveau, Leid een uitdrukking af voor de overgengs= 
wearschijnlijkheid w(n'[N). Sehrijf ter afkorting 
Vod = Ae Hoe luidt nu explicict de Master- 
’ ° 
vergelijking voor de waarsohijnlijkheid P(1,t) dat ie 
N atonen geëxciteerd zijn? dd 
Be Bereken de grootheden ze 
Ek 
A(N) = ) (NLH) WON! IN) $ 
’ ll 
B En 
re: t 2 1 E 
Ao) = (N-N)S WON! [N) | 
GE 
en leid daarmee het rechterlid af van de nomente 
vergelijkingen 
4 8 2N. 
Ù 8 at NP(N,t) = at. . . . | 
8 ie Ì Ti 2 | 
E fe) 2 Me E 
N 





$ 

# 
FS 
In 


C‚ Los de vergelijking voor het eerste moment op. 


Berekon de grootte van de fluotuaties 
Eren we D hi 
Nen 


N, voor evenwicht en zo mogelijk ook voor niet-evenwicht, 





